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Abstract.  Recent  work  has  developed  a  new  framework  for 
the  impact  of  vertical  wind  shear  on  the  intensity  evolution 
of  tropical  cyclones.  A  focus  of  this  framework  is  on  the 
frustration  of  the  tropical  cyclone’s  power  machine  by  shear- 
induced,  persistent  downdrafts  that  flush  relatively  cool  and 
dry  (lower  equivalent  potential  temperature,  9e )  air  into  the 
storm’s  inflow  layer.  These  previous  results  have  been  based 
on  idealised  numerical  experiments  for  which  we  have  de¬ 
liberately  chosen  a  simple  set  of  physical  parameterisations. 
Before  efforts  are  undertaken  to  test  the  proposed  framework 
with  real  atmospheric  data,  we  assess  here  the  robustness  of 
our  previous  results  in  a  more  realistic  and  representative  ex¬ 
perimental  setup  by  surveying  and  diagnosing  five  additional 
numerical  experiments.  The  modifications  of  the  experimen¬ 
tal  setup  comprise  the  values  of  the  exchange  coefficients 
of  surface  heat  and  momentum  fluxes,  the  inclusion  of  ex¬ 
periments  with  ice  microphysics,  and  the  consideration  of 
weaker,  but  still  mature  tropical  cyclones. 

In  all  experiments,  the  depression  of  the  inflow  layer  9e 
values  is  significant  and  all  tropical  cyclones  exhibit  the  same 
general  structural  changes  when  interacting  with  the  imposed 
vertical  wind  shear.  Tropical  cyclones  in  which  strong  down- 
drafts  occur  more  frequently  exhibit  a  more  pronounced  de¬ 
pression  of  inflow  layer  9e  outside  of  the  eyewall  in  our  ex¬ 
periments.  The  magnitude  of  the  9t  depression  underneath 
the  eyewall  early  after  shear  is  imposed  in  our  experiments 
correlates  well  with  the  magnitude  of  the  ensuing  weakening 
of  the  respective  tropical  cyclone.  Based  on  the  evidence  pre¬ 
sented,  it  is  concluded  that  the  newly  proposed  framework  is 
a  robust  description  of  intensity  modification  in  our  suite  of 
experiments. 


1  Introduction 

1.1  A  new  framework  for  intensity  modification  of 
tropical  cyclones  in  vertical  wind  shear 

The  intensity  evolution  of  tropical  cyclones  (TCs)  embedded 
in  an  environmental  flow  with  vertical  wind  shear  poses  an 
important  forecast  problem.  To  date,  the  physical  processes 
that  govern  this  intensity  evolution  are  not  well  understood. 
In  a  recent  publication  (Riemer  et  al.,  2010,  referred  to  as 
RMN10  hereafter)  a  new  framework  for  the  intensity  modi¬ 
fication  of  TCs  in  vertical  wind  shear  was  proposed.  While 
previous  hypotheses  have  focused  almost  exclusively  on  pro¬ 
cesses  above  the  boundary  layer,  RMN10  showed  that  ver¬ 
tical  wind  shear  has  a  significant  impact  on  the  thermody¬ 
namic  properties  of  the  TC’s  inflow  layer1.  Strong  and  per¬ 
sistent,  shear-induced  downdrafts  flush  the  inflow  layer  with 
low  moist-entropy  (or  equivalent  potential  temperature,  9e) 
air  -  “anti-fuel”  for  the  TC  power  machine.  The  replenish¬ 
ment  of  this  air  by  surface  heat  fluxes  is  not  complete  while 
the  air  spirals  towards  the  storm  centre  leading  to  a  reduction 
of  9e  values  within  the  inner-core  updrafts  in  the  azimuthal 
mean.  When  viewed  from  the  simplified  perspective  of  an 
idealised  Carnot  cycle,  a  decrease  in  the  potential  intensity 
of  the  TC  can  be  expected  (Tang  and  Emanuel,  2010). 


'in  the  current  study,  the  term  “inflow  layer”  will  refer  to  the 
layer  of  strong  inflow  associated  with  surface  friction  in  the  lowest 
1-1.5  km  where  the  departure  from  gradient  wind  balance  in  the 
radial  direction  is  generally  a  maximum  (e.g.,  Bui  et  al.  (2009),  their 
Figs.  5  and  6  and  accompanying  discussion). 
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The  downdrafts  that  flush  the  inflow  layer  with  low-0e 
air  are  associated  with  a  quasi-stationary,  azimuthal  wave- 
number  1,  convective  asymmetry  outside  of  the  eyewall.  This 
convective  asymmetry  is  reminiscent  of  the  “stationary  band 
complex”  (Willoughby  et  al.,  1984)  and  for  this  reason  will 
be  referred  to  as  the  SBC  hereafter.  The  SBC  and  the  asso¬ 
ciated  downdraft  pattern  extends  outwards  to  approximately 
150-180  km.  On  account  of  the  swirling  winds,  the  updrafts 
within  the  SBC  rise  along  a  helical  path.  A  simple  kinematic 
model  (Riemer  and  Montgomery,  201 1)  suggests  that,  at  low 
levels  above  the  inflow  layer,  environmental  low-0e  air  ap¬ 
proaches  the  storm  closely  underneath  the  helical  updrafts  in 
the  downshear  and  downshear-left  quadrants2.  Evaporative 
cooling  of  precipitation  falling  out  of  the  SBC,  and  into  this 
unsaturated  air,  most  likely  drives  the  strong  and  persistent, 
vortex-scale  downdrafts  that  cool  and  dry  the  TC’s  inflow 
layer. 

In  RMN10,  evidence  was  presented  that  the  formation  of 
the  SBC  is  connected  to  the  tilt3  of  the  TC  vortex.  When 
subject  to  vertical  wind  shear,  the  tilt  evolution  of  the  TC  vor¬ 
tex  was  shown  to  be  governed  to  zero  order  by  asymmetric- 
balance  dynamics.  Vortex  tilt  can  thus  be  described  in  terms 
of  a  standing  vortex-Rossby-wave  (VRW)  pattern  with  az¬ 
imuthal  and  vertical  wave-number  1  structure.  Even  though 
the  tilt  of  the  inner  core  is  only  10-20  km,  the  VRW  pat¬ 
tern  associated  with  the  tilt  of  the  outer  vortex  extends  out¬ 
wards  to  radii  of  150-200  km  (Fig.  15a  in  RMN10)4.  Af¬ 
ter  vertical  wind  shear  is  suddenly  imposed,  the  TC  settles 
quickly  (approx,  within  3  h)  into  a  quasi-equilibrium  tilt  di¬ 
rection  that  is  to  the  downshear-left.  Thus,  the  positive  vor- 
ticity  anomaly  associated  with  the  standing  wave-number  1 
VRW  pattern  is  found  downshear-left  at  upper  levels  and 
downshear-right  at  low  levels  (Fig.  16  in  RMN10).  Frictional 
convergence  associated  with  this  positive  low-level  vorticity 
anomaly  then  contributes  to  the  forcing  of  the  asymmetric 
convection  within  the  SBC.  The  foregoing  results  suggest  an 
important  connection  between  the  asymmetric -balance  dy¬ 
namics  governing  the  tilt  evolution  of  the  vortex  subject  to 
vertical  shear  and  the  thermodynamic  impact  on  the  inflow 
layer. 

The  findings  of  RMN10  were  derived  from  idealised  nu¬ 
merical  experiments  following  the  setup  of  Bender  (1997) 
and  Frank  and  Ritchie  (2001).  In  this  experimental  setup, 
a  uni-directional  vertical  wind  shear  profile  is  imposed 
suddenly  on  a  mature  model  TC.  In  an  attempt  to  min¬ 


2In  Riemer  and  Montgomery  (2011),  the  distribution  of  9e  was 
analysed  at  2  km  height. 

3  In  RMN10,  as  well  as  in  the  current  study,  the  tilt  is  defined  as 
the  vector  difference  between  the  location  of  the  vorticity  centroids 
at  10  km  and  1  km  height. 

4  A  considerable  tilt  of  the  outer  parts  of  TC-like  vortices  with 

small  tilt  of  the  inner  core  has  been  noted  previously  by  other  au¬ 
thors  also  (see  e.g.  Fig.  6a  in  Jones,  1995;  Fig.  10  in  Reasor  et  al., 
2004). 


imise  complexity  and  to  understand  the  fundamental  pro¬ 
cesses,  RMN10  deliberately  chose  a  simple  set  of  model 
parameterisations.  As  is  the  nature  of  results  derived  from 
such  an  idealised  approach,  uncertainties  exist  whether  the 
newly  proposed  pathway  operates  in  the  real  atmosphere.  Al¬ 
though  some  preliminary  observational  support  was  offered 
in  RMN10  a  more  comprehensive  test  of  the  new  hypothesis 
with  atmospheric  data  awaits  future  studies.  Before  undertak¬ 
ing  such  efforts  we  survey  and  diagnose  here  five  additional 
numerical  experiments  with  some  modifications  of  the  exper¬ 
imental  setup  to  assess  the  robustness  of  RMNlO’s  results. 

1.2  Purpose  of  additional  experiments 

The  experimental  setup  in  RMN10  features  a  simplified 
cloud  microphysics  scheme,  a  likely  overestimation  of  the 
surface  exchange  coefficients  of  momentum  and  enthalpy  at 
high  wind  speed,  and  a  very  high  TC  intensity  at  the  time 
when  shear  is  imposed,  representative  for  a  minority  of  TCs 
in  the  real  atmosphere  only.  The  particular  relevance  of  these 
features  for  RMNlO’s  framework  is  discussed  in  more  detail 
below.  One  goal  of  this  study  is  to  assess  the  robustness  of 
RMNlO’s  results  in  a  more  realistic  and  representative  ex¬ 
perimental  setup. 

Several  environmental  factors  likely  play  a  role  for  the 
evolution  of  TCs  in  vertical  wind  shear  also.  Besides  the 
obvious  importance  of  the  shear  magnitude,  such  factors  in¬ 
clude  the  vertical  profiles  of  environmental  wind  speed  and 
direction  (Zeng  et  al.,  2010;  Wang,  2012),  and  the  environ¬ 
mental  moisture  and  temperature  profiles  (cf.  discussion  in 
Riemer  and  Montgomery  (2011)).  Careful  examination  of 
the  importance  of  these  environmental  profiles  is  beyond  the 
scope  of  this  study  but  constitutes  an  important  topic  for  fu¬ 
ture  research5. 

The  numerical  experiments  in  RMN10  employed  a  micro¬ 
physics  scheme  that  considers  “warm  cloud”  processes  asso¬ 
ciated  with  the  collision-coalescence  process  only  (hereafter 
termed  “warm  rain”).  The  current  study  considers  experi¬ 
ments  with  ice  microphysics  also.  Adding  ice  processes  in¬ 
troduces  the  latent  heat  of  fusion/  sublimation,  new  classes  of 
hydrometeors  with  various  fall  speeds,  and  additional  path¬ 
ways  to  form  precipitation.  Details  of  the  ice  microphysics 
scheme  employed  in  this  study  are  given  in  Sect.  2.2.  Ar¬ 
guably,  the  most  significant  of  the  above  processes  is  the  ad¬ 
ditional  conversion  of  latent  heat  associated  with  the  freezing 
and  melting  of  ice.  Ice  microphysics  produce  stronger  up¬ 
drafts  than  a  warm-rain  scheme  due  to  the  additional  latent 
heat  release  associated  with  freezing.  In  addition,  the  inclu¬ 
sion  of  precipitation  in  the  form  of  ice  may  significantly  en¬ 
hance  the  formation  and  intensity  of  downdrafts  (Srivastava, 
1987). 


5  A  brief  discussion  of  the  potential  consequences  of  the  envi¬ 
ronmental  wind  profile  can  be  found  in  the  authors’  response  to  the 
anonymous  reviewer  (item  1  and  2)  on  the  ACPD  webpage. 
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Previous  studies  have  shown  that  the  inclusion  of  ice  mi¬ 
crophysics  in  numerical  models  may  alter  the  radial  profile  of 
the  swirling  winds  of  TCs  as  compared  to  a  warm-rain  exper¬ 
iment  (e.g..  Lord  et  al.,  1984;  Fovell  et  al.,  2009).  In  general, 
it  is  not  well  understood  how  these  differences  arise6.  Lord 
et  al.  have  noted,  albeit  in  an  axisymmetric  model  frame¬ 
work,  that  their  experiment  with  ice  microphysics  exhibits 
more  mesoscale  convective  features  than  their  experiment 
with  warm  rain.  Fovell  et  al.  have  proposed  a  mechanism  for 
the  change  in  the  radial  profile  of  the  swirling  winds  that  is 
based  on  the  sedimentation  of  ice  particles  from  the  outflow 
anvil,  the  subsequent  sublimentation  of  these  ice  particles, 
and  radiative  feedback.  In  our  suite  of  experiments,  however, 
radiative  processes  are  neglected.  Nevertheless,  potentially 
significant  differences  in  the  TC-vortex  structure  between  the 
warm-rain  and  the  ice  microphysics  experiments  do  arise,  in¬ 
dicating  that  the  mechanism  suggested  by  Fovell  et  al.  is  not 
of  primary  importance.  The  potentially  important  differences 
in  the  TC-vortex  structure  and  their  implication  for  the  in¬ 
teraction  with  vertical  wind  shear  will  be  discussed  in  some 
detail  in  Sect.  3.1. 

Recent  studies  have  highlighted  the  importance  of  the 
representation  of  the  planetary  boundary  layer  in  numeri¬ 
cal  models  for  vortex  intensification  and  structure  changes 
of  simulated  storms  (Braun  and  Tao,  2000;  Smith  and  Thom¬ 
sen,  2010;  Montgomery  et  al.,  2010;  Smith  et  al.,  2012).  The 
exchange  coefficients  for  enthalpy,  Ck,  and  momentum,  Cd, 
play  a  most  prominent  role.  They  govern  the  energy  extrac¬ 
tion  from  the  ocean  surface  and  the  dissipation  of  kinetic 
energy  in  the  frictional  boundary  layer,  respectively.  In  the 
framework  of  RMN10,  the  values  of  Ck  and  Cd  can  be  ex¬ 
pected  to  play  an  important  role  in  the  replenishment  of  the 
depressed  boundary  layer  6t  air  while  spiralling  towards  the 
inner-core  updrafts.  The  replenishment  rate  is  proportional7 
to  Ck  -  The  strength  of  the  inflow  within  the  frictional  bound¬ 
ary  layer  is  determined  in  part  by  Cd-  Smaller  values  of  Co 
may  lead  to  less  radial  inflow  and  thus  to  an  increase  of  the 
time  that  the  air  parcels  spend  under  the  influence  of  sur¬ 
face  fluxes  while  spiralling  inwards.  Assuming  that  Ck  is 
not  changed,  and  that  the  wind  dependence  of  the  fluxes  is 
dominated  by  the  tangential  wind  speed,  the  replenishment 
of  depressed  inflow  layer  6e  values  is  then  more  complete 
than  for  larger  values  of  Co-  Based  on  these  considerations, 
one  would  expect  that  the  shear-induced  flushing  of  the  in¬ 
flow  layer  with  low-0e  air  has  a  more  pronounced  impact  on 
TC  intensity  for  smaller  values  of  Ck  and  larger  values  of 


6In  his  review  of  this  study,  Roger  Smith  has  proposed  that  a 
first-order  explanation  of  how  these  differences  arise  could  be  given 
in  terms  of  the  conventional  (balanced)  model  of  vortex  spin-up, 
noting  that  the  inclusion  of  ice  microphysics  substantially  modifies 
the  diabatic  heating  rate.  We  are  not  aware,  however,  of  published 
studies  that  elaborate  on  this  idea. 

7The  replenishment  rate  is  proportional  to  the  air-sea  thermody¬ 

namic  disequilibrium  also. 


Cd-  The  interpretation  of  the  importance  of  Ck  and  Cd  in 
an  axisymmetric  theory  for  a  TC  in  vertical  wind  shear  by 
Tang  and  Emanuel  (2010)  is  consistent  with  the  foregoing 
considerations. 

Considerable  uncertainties  exist  about  the  values  of  Ck 
and  Cd  at  the  high  wind  speeds  occurring  in  TCs.  Dea¬ 
con’s  formula  for  Cd  was  employed  in  RMN10  and  the  ra¬ 
tio  of  Ck/Cd  was  set  to  unity.  Recent  results  derived  from 
observational  data  (Black  et  al.,  2007;  Zhang  et  al.,  2008) 
strongly  indicate  that  Deacon’s  formula  overestimates  Co  at 
high  wind  speeds  and  that  Ck/Cd  is  considerably  less  than 
unity.  The  observational  data  exhibit  considerable  scatter  and 
are  limited  to  wind  speeds  of  marginal  hurricane  strength8. 
To  address  the  likely  overestimation  of  Ck  and  Ck/Cd  in 
RMN10,  we  employ  in  this  study  values  of  Ck  and  Ck/Cd 
that  are  consistent  with,  but  rather  on  the  lower  end  of  the 
range  of  values  observed  by  Black  et  al.  and  Zhang  et  al. 
Further  details  are  provided  in  Sect.  2.3. 

The  TCs  considered  in  RMN10  were  intense  and  rapidly 
intensifying  at  the  time  when  the  vertical  wind  shear  pro¬ 
file  was  imposed.  In  the  current  study  we  will  consider,  in¬ 
ter  alia,  weaker,  but  still  mature  TCs.  Finally,  our  additional 
experiments  now  exhibit  various  rates  of  intensification  just 
before  vertical  wind  shear  is  imposed.  We  did  not,  however, 
find  a  sensitivity  of  the  subsequent  TC  evolution  in  verti¬ 
cal  shear  to  this  preceding  intensification  rate  in  our  exper¬ 
iments.  Thus,  the  various  rates  of  intensification  just  before 
vertical  wind  shear  is  imposed  will  not  be  considered  here¬ 
after. 

The  next  section  provides  an  overview  of  the  numerical 
model  employed  and  the  additional  experiments  performed, 
an  expanded  description  of  the  ice  microphysics  scheme,  and 
the  values  for  the  surface  exchange  coefficients.  In  Sect.  3 
we  compare  the  evolution  of  the  quintessential  ingredients 
of  RMNlO’s  newly  proposed  framework  in  the  different  nu¬ 
merical  experiments.  One  experiment  with  ice  microphysics 
is  examined  in  further  detail  in  Sect.  4,  including  an  estimate 
of  the  spindown  timescale  of  the  simulated  TC-vortex  based 
on  simple  axisymmetric  theory.  Section  5  contains  our  con¬ 
clusions. 

2  Model  setup  and  additional  experiments 
2.1  A  brief  model  description 

For  the  numerical  experiments  we  employ  the  Regional 
Atmospheric  Modeling  System  (RAMS),  developed  at 
Colorado  State  University  (Pielke  et  al.,  1992;  Cotton 
et  al.,  2003).  The  RAMS  is  a  three-dimensional,  non¬ 
hydrostatic  numerical  modeling  system  comprising  time- 
dependent  equations  for  velocity,  non-dimensional  pressure 

8  Some  significant  strides  towards  determining  these  coefficients 
at  major  hurricane  wind  conditions  have  been  made  recently  by  Bell 
et  al.  (2012). 


www.atmos-chem-phys.net/13/327/2013/ 


Atmos.  Chem.  Phys.,  13,  327-346,  2013 


330 


M.  Riemer  et  al.:  Thermodynamic  impact  of  vertical  wind  shear  on  the  inflow  layer 


perturbation,  ice-liquid  water  potential  temperature  (Tripoli 
and  Cotton,  1981),  and  mixing  ratios  of  total  water  and  hy¬ 
drometeors. 

The  experiments  are  carried  out  on  an  f-plane  at  15°  N  and 
with  a  horizontal  grid  spacing  of  5  km.  The  model  does  not 
include  ocean  feedback.  The  time-invariant  sea  surface  tem¬ 
perature  is  set  to  28.5°  C.  A  standard  first-order  turbulence 
scheme  developed  by  Smagorinsky  (1963)  is  used  for  sub- 
grid-scale  mixing,  with  modifications  by  Lilly  (1962)  and 
Hill  (1974)  that  enhance  diffusion  in  unstable  conditions  and 
reduce  diffusion  in  stable  conditions.  For  simplicity,  radiative 
processes  are  neglected.  The  bulk  aerodynamic  formulas  are 
used  to  calculate  momentum  fluxes  and  sensible  and  latent 
heat  fluxes  from  the  sea  surface.  The  values  of  the  exchange 
coefficients  are  given  in  Sect.  2.3.  We  employ  either  a  warm- 
rain  microphysics  scheme  based  on  that  of  Kessler  (1969) 
or  the  RAMS  standard  ice  microphysics  scheme  (Sect.  2.2). 
The  mixing  ratio  equations  for  all  prognosed  hydrometeors 
include  advection  by  the  three-dimensional  flow.  A  more 
comprehensive  description  of  RAMS,  the  surface  interac¬ 
tion  parameterisation,  and  the  warm-rain  scheme  is  given  in 
RMN10. 

2.2  Ice  microphysics  in  RAMS 

The  RAMS  ice  microphysics  scheme  has  eight  water  cate¬ 
gories:  vapour,  cloud  droplets,  rain,  pristine  ice,  snow,  ag¬ 
gregates,  graupel  and  hail  (Walko  et  al.,  1995;  Meyers  et  al., 
1997).  The  liquid  phase  categories,  cloud  droplets  and  rain, 
may  be  supercooled.  Pristine  ice,  snow  and  aggregates  are 
assumed  to  be  completely  frozen,  whereas  graupel  and  hail 
are  mixed  phase  categories,  capable  of  consisting  of  ice  only, 
or  a  mixture  of  ice  and  liquid.  All  hydrometeor  categories  ex¬ 
cept  cloud  droplets  are  assumed  large  enough  to  fall.  Cloud 
droplets  and  pristine  ice  are  the  only  categories  to  nucleate 
from  vapour.  All  other  categories  form  from  existing  hy¬ 
drometeors,  but  once  formed  may  grow  further  by  vapour 
deposition.  Pristine  ice  crystals  form  from  heterogeneous  nu- 
cleation,  homogeneous  nucleation,  and  secondary  produc¬ 
tion  that  involves  splintering  of  riming  crystals.  The  snow 
category  is  defined  in  RAMS  as  relatively  large  ice  crys¬ 
tals,  which  have  grown  from  the  pristine  ice  category  due 
to  vapour  deposition  and  may  continue  to  grow  by  riming. 
Aggregates  are  defined  as  ice  particles  that  have  formed  by 
collision  and  coalescence  of  pristine  ice,  snow  and/or  aggre¬ 
gates.  Snow  and  aggregates  retain  their  identity  with  moder¬ 
ate  amounts  of  riming,  but  convert  to  graupel  if  the  amount  of 
riming  is  large.  In  addition  to  forming  by  heavy  riming,  grau¬ 
pel  can  also  form  by  partial  melting  of  snow  or  aggregates. 
Graupel  is  allowed  to  carry  up  to  half  its  mass  in  liquid.  If 
the  percentage  becomes  larger,  by  either  melting,  riming  or 
collection  by  rain,  a  graupel  particle  is  re-categorised  as  hail. 
Hail  is  formed  by  freezing  of  rain  drops  or  by  riming  or  par¬ 
tial  melting  of  graupel.  Hail  is  allowed  to  be  a  large  fraction 
liquid  and  converts  to  rain  when  completely  liquid. 


The  hydrometeor  size  distribution  in  RAMS  is  represented 
by  a  generalised  Gamma  function.  A  user  specified  shape  pa¬ 
rameter  controls  the  width  of  the  distribution.  Also,  a  scaling 
diameter  is  often  specified  for  a  one-moment  scheme  that  is 
related  to  the  mean  diameter  of  the  distribution,  with  the  par¬ 
ticular  relation  depending  on  the  chosen  spectral  parameters. 
For  this  study  a  one-moment  scheme  is  employed  for  all  cate¬ 
gories  except  cloud  droplets  and  pristine  ice.  A  scaling  diam¬ 
eter  is  specified  and  the  number  concentration  is  diagnosed 
using  the  predicted  mixing  ratio  of  the  hydrometeor  scheme. 

A  heat  budget  equation  is  formulated  for  rain,  graupel,  and 
hail  hydrometeors.  This  allows  the  temperature  of  a  hydrom¬ 
eteor  to  differ  from  that  of  the  surrounding  ambient  air  tem¬ 
perature  due  to  latent  heat  release  or  absorption  in  the  hy¬ 
drometeor,  and  sensible  heating  by  collisions  with  other  hy¬ 
drometeors.  The  temperature  of  a  hydro  meteor  often  differs 
substantially  from  that  of  the  air,  which  can  significantly  in¬ 
fluence  the  rates  of  heat  and  vapour  diffusion,  as  well  as  the 
amount  of  sensible  heat  transfer  that  occurs  in  coalescence 
of  hydro  meteors.  Collision  and  coalescence  of  hydrometeors 
is  governed  by  the  stochastic  equation  derived  by  Verlinde 
et  al.  (1990). 

Gravitational  settling  of  hydrometeors  causes  them  to  fall 
relative  to  air.  Terminal  velocities  increase  as  a  power  of  the 
droplet  diameter,  and  are  based  on  empirical  data.  Pristine 
ice  has  very  small  fall  speeds.  Snow  crystals  fall  faster  than 
pristine  ice,  but  still  relatively  slowly.  Aggregates  are  large 
crystals  that  have  moderate  fall  speeds.  Graupel  and  hail  have 
high  fall  speeds.  Sedimentation  in  RAMS  uses  a  Lagrangian 
scheme  to  transport  the  mixing  ratio  from  any  given  grid  cell 
to  a  lower  height  in  the  vertical  column,  so  that  a  proportion 
is  transferred  to  the  lower  adjacent  grid  cell. 

2.3  Modification  of  surface  exchange  coefficients 

Deacon’s  formula  for  the  dependence  of  the  drag  exchange 
coefficient  C d  on  wind  speed  is  given  by 

CD  =  1.1  x  10“3+4x  10“5|V|i,  (1) 

with  |  V 1 1  being  the  wind  speed  in  m  s- 1  at  the  lowest  model 
level  (here:  49  m).  In  RMN10  we  used  Deacon’s  formula  to¬ 
gether  with  C'k  —  Cd-  For  reference,  this  formula  furnishes 
values  of  the  exchange  coefficients  of  1.5  x  10-3,  2.1  x  10-3, 
and  3.1  x  10-3  for  10ms_1,  25ms-1,  and  50ms-1,  respec¬ 
tively. 

Zhang  et  al.  (2008)  found  that  “77ie  exchange  coefficient  of 
enthalpy  flux  shows  no  significant  dependence  on  wind  speed 
up  to  hurricane  force  with  a  value  of  1 .16  x  70-3.”  and  that 
“77ze  average  ratio  of  the  Ck/Co  values  is  0.63,  ...”.  We 
therefore  use  these  values  and  set 

Ck  =  1.16  x  10-3  =  constant  (2) 

and 

CK/ Cd  —  0.63.  (3) 
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Table  1.  Overview  of  the  experiments  that  are  considered  in  this  study.  The  experiments  RMNgg  and  jgRMNgjj  are  described  in  detail  in 
RMN10.  In  column  Ck,  “Deacon”  denotes  that  Deacon’s  formula  (see  Sect.  2.3)  is  used  to  determine  the  exchange  coefficients.  The  last 
column  gives  the  TC's  intensification  rate  just  before  vertical  shear  is  imposed. 


Experiment 

“initial” 

intensity  [ms-1] 

CK 

Ck/Cd 

microphysics 

shear 

[ms-1 /12  km] 

preceding  int. 
rate  [m  s-1/h] 

rmn68 

68 

Deacon 

1 

warm 

15 

1.8 

10RMN68 

68 

Deacon 

1 

warm 

10 

1.8 

RMN54 

54 

Deacon 

1 

warm 

15 

0.0 

CBLAST68 

68 

1.16  x  10-3 

0.65 

warm 

15 

1.1 

cblast54 

54 

1.16  x  10-3 

0.65 

warm 

15 

0.0 

ice68 

68 

Deacon 

1 

ice 

15 

2.5 

1qICE68 

68 

Deacon 

1 

ice 

10 

2.5 

Experiments  that  employ  these  values  are  referred  to  as 
CBLAST  below.  Note  that  Deacon’s  formula  and  Ck/Cd  — 
1  furnishes  values  of  Ck  and  Cd  that,  for  |  V|i  =  50 ms-1, 
are  approx.  2.7  and  4.2  times  higher,  respectively,  than  the 
corresponding  CBLAST  values. 

2.4  Summary  of  experiments 

A  summary  of  the  experiments  examined  in  this  study  is 
given  in  Table  1.  We  refer  to  experiments  that  employ 
the  same  parameterisations  as  in  RMN10,  i.e.  warm  rain 
and  Deacon’s  formula  for  the  exchange  coefficients  (with 
Ck/Cd  —  1).  as  RMN.  The  CBLAST  experiments  employ 
warm  rain  also.  Experiments  that  include  ice  microphysics 
are  referred  to  as  ICE.  ICE  experiments  use  Deacon’s  for¬ 
mula  with  Ck/Cd  —  L  Besides  the  modifications  that  are 
mentioned  explicitly,  the  experimental  setup  is  the  same  as 
in  RMN  10. 

The  shear  profile  employed  in  this  study  has  a  sine  struc¬ 
ture  in  the  vertical  with  zero  winds  at  the  surface  and  max¬ 
imum  easterly  winds  at  12  km  and  above.  We  will  follow 
the  convention  of  the  TC  community  to  define  vertical  wind 
shear  as  the  vector  difference  between  the  winds  at  upper 
and  lower  levels.  Here,  we  use  the  vector  difference  between 
12  km  and  the  surface.  The  shear  magnitude  is  thus  given 
in  ms-1.  Unless  otherwise  noted,  we  will  consider  TCs  in 
vertical  wind  shear  with  a  magnitude  of  15  ms-1.  The  shear 
profile  is  imposed  at  a  high  intensity9  of  68  ms-1,  after  the 
TCs  have  spun  up  in  quiescent  environment.  These  experi¬ 
ments  are  referred  to  as  RMNg8,  CBLAST^s,  and  ICEgg,  re¬ 
spectively.  The  intensity  at  the  time  when  shear  is  imposed 
will  be  referred  to  as  “initial  intensity”  below.  The  conven¬ 
tions  used  in  this  study  follow  those  in  RMN10.  Experiment 
RMN68  is  the  “15mps”  case  of  RMN10  and  constitutes  our 
reference  experiment. 

We  have  performed  experiments  with  a  weaker  ini¬ 
tial  intensity  of  54 ms-1  also,  referred  to  as  RMN54  and 


9  Our  storm  intensity  metric  is  the  maximum  azimuthally- 
averaged  tangential  wind  speed  at  1  km  height. 


CBLAST54.  These  experiments  use  the  same  spinup  runs  as 
the  respective  experiments  with  initial  intensity  of  68  ms-1, 
but  with  vertical  shear  imposed  earlier  during  the  inten¬ 
sification.  Furthermore,  an  ICE  experiment  has  been  per¬ 
formed  with  10ms_1  vertical  shear  (ipICEgg)  which  will 
be  compared  briefly  to  the  “lOmps”  case  of  RMN10  (here: 
ioRMN6g). 

3  Assessing  the  robustness  of  the  new  framework 

We  shall  focus  on  the  key  features  of  RMNlO’s  new  frame¬ 
work  for  intensity  modification  in  vertical  wind  shear:  the 
vortex  tilt,  the  occurrence  and  location  of  the  SBC,  the  flush¬ 
ing  of  the  inflow  layer  with  low-0e  air,  the  ensuing  depres¬ 
sion  of  inflow  layer  Be,  and  the  associated  weakening  of  the 
TC.  The  general  characteristics  of  these  key  features  will  be 
found  to  be  common  in  all  experiments. 

As  can  be  expected,  differences  arise  in  the  specific  char¬ 
acteristics  of  these  key  features.  For  some  of  these  differ¬ 
ences  it  is  very  difficult  to  establish  the  causal  link  to  the 
changes  in  the  experimental  setup.  Most  prominently,  differ¬ 
ences  between  the  evolution  of  the  TCs  in  RMNgg  and  ICEf,g 
are  partly  due  to  the  enhancement  of  downdrafts  by  ice  mi¬ 
crophysics  (Sect.  3.7),  but  very  likely  due  also  to  changes  in 
the  radial  profile  of  the  swirling  winds  at  the  time  when  verti¬ 
cal  wind  shear  is  imposed.  Disentangling  the  individual  con¬ 
tributions  requires  a  more  complete  understanding  of  the  pro¬ 
cesses  at  work,  and  extended  numerical  experimentation  and 
diagnostic  analyses,  preferably  in  an  ensemble-based  frame¬ 
work.  A  more  comprehensive  examination  of  such  differ¬ 
ences  is  thus  deferred  to  future  research. 

3.1  Importance  of  the  radial  profile  of  the  swirling 
winds 

3.1.1  Impact  on  the  resiliency  of  TC-like  vortices 

It  is  generally  believed  that  the  presence  of  vertical  wind 
shear  increases  the  interaction  between  a  TC  and  dry  environ¬ 
mental  air  (Riemer  and  Montgomery,  2011,  and  references 
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therein).  Based  on  a  simple  kinematic  model,  Riemer  and 
Montgomery  (201 1)  argue  that  a  TC  with  a  broad  radial  pro¬ 
file  of  the  swirling  winds  should  have  a  stronger  ability  to 
isolate  itself  from  the  detrimental  interaction  with  dry  envi¬ 
ronmental  air  than  a  smaller-scale  TC.  Furthermore,  the  ra¬ 
dial  profile  of  the  swirling  winds  has  a  profound  impact  on 
the  dynamic  resiliency  of  TC-like  vortices,  i.e.  the  ability  to 
maintain  vertical  coherence  in  the  presence  of  differential  ad- 
vection  by  the  vertical  shear  flow  (e.g.,  Jones,  1995;  Reasor 
et  al.,  2004).  These  studies  investigated  the  resiliency  of  ini¬ 
tially  barotropic,  dry  vortices  and  are  reviewed  briefly  below. 

Jones  (1995)  examined  vortex  evolution  in  terms  of  the 
formation  and  subsequent  interaction  of  potential  vorticity 
(PV)  anomalies.  PV  anomalies  were  defined  as  deviations 
from  the  azimuthal  mean.  Due  to  the  differential  advection 
by  the  vertical  shear  flow,  PV  anomalies  representing  the  tilt 
of  the  vortex  form  at  lower  and  upper  levels.  These  anoma¬ 
lies  mutually  interact  through  the  vertical  penetration  of  their 
“induced”  quasi-horizontal  circulations,  causing  the  vortex 
to  precess.  This  precession  helps  the  vortex  to  counteract  the 
differential  advection  by  the  vertical  wind  shear.  The  verti¬ 
cal  interaction  of  these  PV  anomalies  depends  on  the  Rossby 
penetration  depth.  The  penetration  depth  increases  with  the 
strength  of  the  background  rotation  and  with  the  horizontal 
scale  of  the  interacting  anomalies10.  Jones  noted  that  it  is 
not  entirely  clear  how  to  quantify  the  penetration  depth  for 
a  tilted  TC-like  vortex  but  concluded  that  intense  and  broad 
vortices  are  more  resilient  than  their  weaker,  smaller-scale 
counterparts. 

Reasor  et  al.  (2004)  re-visited  the  problem  of  vortex  re¬ 
siliency  and  re-alignment  using  VRW  concepts  and  models. 
Their  results  are  in  qualitative  agreement  with  Jones’s  find¬ 
ing  that  vortex  (and  environment)  parameters  that  increase 
the  penetration  depth  improve  the  resistence  of  a  vortex  to 
tilting  by  an  imposed  vertical  shear  flow.  In  addition,  Reasor 
et  al.  demonstrated  that  vortex  precession  is  not  the  only  ele¬ 
ment  in  the  resiliency  of  TC-like  vortices.  An  inviscid  damp¬ 
ing  mechanism  of  the  vortex-tilt  mode  exists,  associated  with 
PV  mixing  at  the  critical  radius  of  the  near-discrete  VRW 
that  represents  the  vortex  tilt.  The  critical  radius  is  typically 
located  in  the  “skirt”  region  of  the  vortex,  i.e.  in  the  region 
of  small  PV  values  outside  of  the  high  PV  values  in  the  in¬ 
ner  core  (approx,  outside  50-70 km  in  Fig.  la).  The  invis¬ 
cid  damping  mechanism  depends  on  the  radial  gradient  of 
the  skirt  PV  at  the  critical  radius  (Schecter  and  Montgomery, 
2003).  Inviscid  damping,  and  thus  re-alignment  of  the  vortex, 
occurs  when  the  radial  PV  gradient  in  this  region  is  negative. 
The  damping  rate  increases  with  the  magnitude  of  the  gra¬ 
dient.  For  a  positive  skirt-PV  gradient,  a  tilt  instability  may 
occur.  The  results  of  Reasor  et  al.  clearly  demonstrate  that 
vortex  resiliency  may  crucially  depend  on  the  radial  profile 
of  the  swirling  winds. 


10The  penetration  depth  increases  also  with  a  decrease  in  static 
stability. 


Fig.  1.  (a)  Radial  profiles  of  azimuthal  mean  vorticity  at  2  km  (thin) 
and  azimuthal  mean  tangential  wind  at  the  radius  of  maximum 
winds  at  1  km  height  (thick)  at  the  time  when  vertical  wind  shear  is 
imposed  in  CBLAST54  (red),  CBLASTgg  (blue),  and  ICEgg  (grey), 
respectively,  (b)  Radial  vorticity  profiles  as  in  (a)  but  zoomed  in  on 
the  “skirt”  region  between  50  km  and  100  km. 


The  radial  profile  of  the  swirling  winds  has  been  confirmed 
to  play  a  role  in  the  intensity  evolution  of  vertically  sheared 
TCs  in  idealised,  moist  numerical  experiments  (Wong  and 
Chan,  2004).  In  these  experiments  the  weakening  of  smaller- 
scale  TCs  is  found  to  be  more  pronounced  than  that  of  their 
broader-scale  counterparts.  Wong  and  Chan  based  their  inter¬ 
pretation  of  this  difference  in  the  intensity  evolution  on  the 
results  of  Jones  (1995)  derived  from  dry  TC-like  vortices.  To 
us,  however,  it  is  not  clear  what  the  underlying  processes  are 
that  govern  the  intensity  evolution  in  Wong  and  Chan’s  ex¬ 
periments.  In  particular,  the  role  of  moist  processes  and/or 
the  associated  secondary  circulation  remains  unclear. 

In  the  real  atmosphere,  an  indication  of  a  higher  suscepti¬ 
bility  of  smaller-scale  TCs  to  vertical  shear  has  been  found 
by  DeMaria  (1996)  based  on  a  regression  analysis. 

3.1.2  Diagnosed  radial  vortex  structure  in  our 
experiments 

Because  the  radial  profile  of  the  swirling  winds  may  play 
a  prominent  role  in  the  evolution  of  TCs  in  vertical  wind 
shear,  we  begin  the  examination  of  our  experiments  with  the 
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Fig.  2.  (a)  Intensity  time  series  for  the  RMN  no-shear  experiment  (red),  RMN54  (light  red),  RMNgg  (dark  red),  the  CBLAST  no-shear 
experiment  (blue),  CBLAST54  (light  blue),  and  CBLASTgg  (dark  blue).  The  colour-coded  arrows  denote  the  time  when  vertical  wind  shear 
is  imposed  in  the  respective  experiment,  (b)  As  in  (a),  but  for  the  RMN  no-shear  experiment  (red),  RMNgg  (dark  red),  jgRMNgg  (light  red), 
the  ICE  no-shear  experiment  (black),  ICEgg  (dark  grey),  and  ujICEgg  (light  grey),  (c)  Time  series  of  the  difference  between  the  intensity  in 
RMN54  (light  red),  RMNgg  (dark  red),  CBLAST54  (light  blue),  and  CBLASTgg  (dark  blue)  and  the  respective  no-shear  experiment.  Time 
Oh  denotes  the  time  when  vertical  wind  shear  is  imposed  in  the  respective  experiment,  (d)  As  in  (c),  but  for  RMNgg  (dark  red),  igRMNgg 
(red),  ICE6g  (dark  grey),  and  iqICE68  (light  grey). 


comparison  of  the  vortex  structure  at  the  time  when  vertical 
wind  shear  is  imposed.  The  radial,  azimuthal-mean  profiles 
of  the  swirling  winds  at  this  time  are  illustrated  in  Fig.  1.  In 
CBLAST68,  the  radius  of  maximum  winds  (RMW)  is  found 
at  40km  (Fig.  la).  The  vorticity  maximum  is  35  x  10-4  s-1 
at  a  radius  of  20  km.  This  vortex  structure  is  very  similar  as  in 
RMN68  (shown  in  Fig.  4  in  RMN10)11.  For  the  weaker  TC  in 
CBLAST54,  the  vorticity  structure  outside  a  radius  of  35  km 
is  almost  identical  to  CBLASTgg.  Inside  this  radius,  vortic¬ 
ity  is  virtually  constant  with  a  value  of  approx.  20  x  1 0  4  s  1 . 
The  corresponding  RMW  is  found  at  48  km.  The  initial  radial 
structure  of  the  TC  in  RMN54  is  very  similar  (not  shown). 

The  vorticity  structure  of  the  TC  in  ICE^g  exhibits  some 
notable  differences  as  compared  to  the  warm-rain  experi¬ 
ments  CBLASTfig  and  RMN'(,g.  The  vorticity  maximum  is 
found  inside  a  radius  of  10  km  with  values  of  55  x  10~4  s  , 
i.e.  approx.  60  %  higher  than  in  CBLAST68-  Radially  outside 
of  this  maximum,  the  high  vorticity  values  decrease  rapidly. 


1 1 A  potential  impact  of  barotropic  instability  associated  with  the 
reversal  of  the  radial  vorticity  gradient  in  the  inner  core  on  the  in¬ 
tensity  evolution  of  the  sheared  TCs  in  our  experiments  is  arguably 
negligible  (see  Sect.  3.1  in  RMN10). 


Corresponding  to  the  smaller  radius  of  maximum  vorticity, 
the  RMW  is  smaller  (located  at  30  km)  than  in  CBLASTgg 
and  RMN68  also. 

As  discussed  above,  the  resiliency  of  the  TC  vortex  may 
crucially  depend  on  the  gradient  of  (potential)  vorticity  at  the 
critical  radius  of  the  (near-discrete)  VRW  mode  that  repre¬ 
sents  the  vortex  tilt.  The  critical  radius  is  found  where  the 
azimuthal  propagation  of  this  VRW,  i.e.  the  vortex  preces¬ 
sion,  equals  the  speed  of  the  swirling  winds.  The  precession 
rate,  and  thus  the  critical  radius,  in  our  experiments  is  diffi¬ 
cult  to  determine.  After  the  vortices  have  re-aligned,  a  short 
precession  period  of  2-3  h  is  indicated  in  all  experiments 
(see  Fig.  4  below12).  Using  the  values  of  the  axisymmetric 
swirling  winds  at  1  km,  as  depicted  in  Fig.  la,  our  rough  es¬ 
timate  of  the  critical  radius  is  then  50-90  km. 

In  ICE,  the  vorticity  gradient  is  near-zero  between  50  km 
and  65  km  while  in  CBLAST68  and  RMN54  the  gradient  is 
negative  and  of  large  magnitude  (Fig.  lb).  Between  65  km 
and  100  km  the  radial  vorticity  profiles  are  broadly  similar. 
At  larger  radii  (not  shown)  the  average  vorticity  gradient  is 


12Note  that  the  temporal  resolution  of  the  data  in  Fig.  4  is  1  h.  It 
is  thus  most  likely  that  the  vortex  precession  is  undersampled. 
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positive  in  RMN54  (between  120  km  and  200  km)  and  ICE 
(between  150-200km).  In  contrast,  in  CBLAST,  vorticity 
decreases  monotonically  out  to  200  km  radius.  It  is  question¬ 
able,  however,  whether  these  differences  at  a  radius  consider¬ 
ably  larger  than  our  estimate  for  the  critical  radius  play  a  role 
for  the  TC-vortex  evolution  in  our  experiments. 

From  the  discussion  in  this  subsection  it  is  clear  that  dif¬ 
ferences  in  the  evolution  of  the  TCs  when  interacting  with 
the  vertical  shear  flow  in  the  various  experiments  potentially 
arise  due  to  differences  in  the  radial  vorticity  profiles  that  ex¬ 
ist  before  shear  is  imposed.  The  weaker  vortices  in  RMN54 
and  CBLAST54,  as  well  as  the  smaller-scale  vortex  in  ICEes, 
can  be  expected  to  be  more  susceptible  to  tilting  than  their 
stronger,  respectively  broader-scale  counterparts  in  RMN68 
and  CBLAST68.  In  rCEf,8-  the  susceptibility  to  tilting  is  po¬ 
tentially  enhanced  by  the  near-zero  gradient  of  the  skirt  vor¬ 
ticity  inside  of  65  km.  In  Sect.  3.3  it  will  be  shown  that  the 
tilt  evolution  is  broadly  consistent  with  these  expectations. 
Furthermore,  based  on  the  results  of  Wong  and  Chan  (2004), 
summarized  at  the  end  of  Sect.  3.1.1,  a  more  pronounced  in¬ 
tensity  decrease  in  ICE68  can  be  expected  due  to  the  smaller 
radial  scale  of  the  vortex  alone.  The  latter  point  needs  to  be 
borne  in  mind  when  attempting  to  attribute  differences  be¬ 
tween  ICE68  and  RMN68  to  their  different  representation  of 
cloud  microphysical  processes. 

3.2  Intensity  evolution 

An  overview  of  the  intensity  evolution  in  our  experiments 
is  presented  in  Fig.  2.  First,  a  brief  account  of  the  intensity 
evolution  in  the  quiescent  environment,  i.e.  without  verti¬ 
cal  shear,  is  given.  The  intensity  evolution  in  RMN  is  de¬ 
scribed  in  some  detail  in  RMN10.  In  CBFAST  (Fig.  2a), 
the  incipient  TC  slowly  intensifies  to  25 ms-1  in  the  first 
20  h.  A  period  of  rapid  intensification  ensues  during  which 
the  TC  reaches  an  intensity  of  53 ms-1  at  33 h.  In  the  sub¬ 
sequent  20  h,  the  TC  consolidates  and  slightly  intensifies  to 
56  ms-1.  Then,  the  TC  continues  to  intensify  at  a  moderate 
rate  until  it  reaches  a  quasi-steady  intensity  of  80-85  ms-1  at 
80 h.  In  ICE  (Fig.  2b),  the  TC  intensifies  to  approx.  40 ms-1 
at  30  h  after  a  short  gestation  period.  Very  rapid  intensifica¬ 
tion  ensues  and  the  intensity  reaches  90 ms-1  at  50 h.  The 
maximum,  quasi-steady  intensity  of  95-100 ms-1  is  reached 
around  72  h. 

In  all  experiments,  intensity  decreases  after  vertical  wind 
shear  is  imposed.  As  in  RMN10,  all  TCs  are  resilient  to 
the  detrimental  impact  of  vertical  wind  shear.  After  a  period 
of  weakening  the  TCs  re-intensify;  particularly  rapidly  in 
RMN5413.  The  re-intensification  process  in  this  experimen¬ 
tal  configuration  is  discussed  in  some  detail  in  RMN10  (their 
Sect.  6). 


13  Because  the  focus  of  this  study  is  on  the  structural  changes  of 
the  TCs  during  the  weakening  phase,  this  experiment  was  discon¬ 
tinued  at  84  h. 


To  facilitate  the  intercomparison  of  the  intensity  evolu¬ 
tion  in  the  various  experiments,  the  intensity  difference  from 
the  respective  no-shear  experiment  is  presented  in  Fig.  2c,  d. 
Time  0  h  in  these  figures  denotes  the  time  when  vertical  wind 
shear  is  imposed  in  the  respective  experiment.  In  all  warm- 
rain  experiments  (Fig.  2c),  the  characteristics  of  this  relative 
weakening  are  similar  in  the  first  20  h.  The  only  exception  be¬ 
ing  that  the  weakening  in  CBFAST54  is  delayed  by  4h.  This 
detail  will  be  examined  in  Sect.  3.6.2.  After  20  h,  the  exper¬ 
iments  with  weaker  TCs  (RMN54  and  CBFAST54)  continue 
to  weaken.  In  contrast,  CBEAST^s  maintains  its  intensity 
difference  for  the  next  20  h  and  RMN68  re-intensifies.  Based 
on  these  plots  alone  it  might  therefore  be  concluded  that  the 
weaker  storms  are  more  susceptible  to  shear.  However,  this 
is  not  necessarily  the  case  here.  Note  that  20  h  after  shear  is 
imposed  in  RMN54  and  CBFAST54  the  respective  reference 
TCs  in  a  quiescent  environment  are  still  rapidly  intensifying. 
In  contrast,  the  reference  TCs  in  RMN  6  8  and  CBLAS’IViS  in¬ 
tensify  only  slowly  20  h  after  vertical  wind  shear  is  imposed. 
Thus,  the  differences  in  the  relative  intensity  evolution  partly 
arise  due  to  the  differences  in  the  contemporaneous  evolution 
of  the  reference  TC.  This  brief  discussion  illustrates  one  of 
the  inevitable  challenges  to  interpreting  different  TC  simula¬ 
tions  in  vertical  wind  shear. 

The  comparison  of  the  RMN  and  ICE  shear  experiments 
is  less  complicated  because  the  intensity  evolution  of  the 
reference  runs  in  quiescent  environment  is  similar  (Fig.  2b, 
d).  A  comparison  of  the  intensity  evolution  in  loICEgs  and 
ICE68,  respectively,  clearly  shows  that  the  TCs  in  the  ICE 
experiments  are  more  susceptible  to  shear  than  their  counter¬ 
parts  in  the  warm-rain  experiments.  The  differences  between 
the  relative  weakening  in  loICEgs  and  ioRMN^s,  i.e.  the 
experiments  with  imposed  vertical  shear  of  10ms_1,  is  on 
average  approx.  10  ms-1  before  16  h,  and  on  average  ap¬ 
prox.  5  ms-1  between  16  h  and  36  h  (Fig.  2d).  After  this 
time,  it  is  no  longer  evident  that  the  weakening  in  joICEes 
is  more  pronounced  than  in  k)RMN68- 

For  15  ms-1  shear,  the  difference  in  the  intensity  evolu¬ 
tion  between  the  ice  and  warm-rain  experiments  is  striking. 
The  maximum  relative  weakening  in  rCE(,s  is  —60  ms-1 
(Fig.  2d).  The  TC  reaches  a  minimum  intensity  of  35  ms-1 
(Fig.  2b).  For  RMN^,  the  respective  values  are  —25ms-1 
and  63ms-1.  Notwithstanding  this  distinct  intensity  evolu¬ 
tion,  it  will  be  shown  below  that  the  structural  changes  in 
ICE(,8  and  RMN^  are  very  similar.  In  particular,  the  pro¬ 
nounced  weakening  in  ICEes  is  associated  with  by  far  the 
most  significant  depression  of  inflow  layer  6e  (Sects.  3.6 
and  4.1).  The  evolution  in  ICEgs  is  thus  consistent  with  the 
results  of  RMN10,  i.e.  the  shear-induced  frustration  of  in¬ 
flow  layer  9e  plays  a  crucial  role  in  the  intensity  evolution  in 
ICE68-  Furthermore,  it  will  be  argued  in  Sect.  4.2  that  this 
rapid  weakening  is  consistent  with  the  frictional  spindown  of 
an  axisymmetric  vortex  due  to  divergence  above  the  inflow 
layer  associated  with  a  significant  reduction  of  the  inner-core 
convective  mass  flux. 
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Fig.  3.  (a)  Time  series  of  the  tilt  magnitude  for  RMN54  (light  red).  RMN^g  (dark  red),  CBLAST54  (light  blue),  and  CBLASTgg  (dark  blue). 
Time  0  h  denotes  the  time  when  vertical  wind  shear  is  imposed  in  the  respective  experiment.  Vortex  tilt  is  defined  as  the  vector  difference 
between  the  vorticity  centroids  at  10  km  and  1  km  height.  For  visual  clarity,  the  3h  running  mean  of  the  time  series  is  shown,  (b)  As  in  (a) 
but  for  RMNgg  (dark  red),  ipRMNgg  (light  red),  ICEgg  (dark  grey),  and  iqlCEgg  (light  grey).  Note  that  the  tilt  magnitude  for  ICEgg  goes 
off  the  scale  between  33  h  and  38  h  reaching  a  maximum  value  of  195  km  (not  shown). 


Fig.  4.  Time  series  of  the  tilt  direction  for  RMNgg  (a),  CBLASTgg  (b),  RMN54  (c),  CBLAST54  (d),  ICEgg  (e),  and  njICEgg  (f).  Time  0  h 
denotes  the  time  when  vertical  wind  shear  is  imposed  in  the  respective  experiment.  Vortex  tilt  is  defined  as  the  vector  difference  between  the 
vorticity  centroids  at  10  km  and  1  km  height.  Dashed  line  denotes  the  3  h  running  mean.  Note  that  the  temporal  resolution  (1  h)  of  the  data  is 
too  coarse  to  resolve  a  putative  high-frequency  precession  of  the  TC  vortex. 


3.3  Evolution  of  vortex  tilt 

Following  RMN10,  the  tilt  of  the  vortex  is  defined  here  as 
the  vector  difference  between  the  vorticity  centroids  at  10  km 
and  1  km  height.  Details  on  the  calculation  of  the  vorticity 
centroids  are  given  in  RMN10.  Here,  the  tilt  magnitude  and 


the  tilt  direction  in  our  experiments  are  depicted  in  Figs.  3 
and  4,  respectively. 

A  detailed  examination  of  the  tilt  evolution  in  RMN(,g 
was  given  in  Sect.  5.1  of  RMN10.  There,  evidence  was  pre¬ 
sented  that  the  tilt  evolution  is  governed  to  zero  order  by 
asymmetric -balance  dynamics.  The  TC  in  RMN68  settles 
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into  a  quasi-equilibrium  tilt  direction  to  the  downshear-left 
(south)  shortly  after  shear  is  imposed  (approx,  within  3  h). 
A  quasi-equilibrium  tilt  direction  may  be  achieved  when 
the  cyclonic  precession  of  the  tilted  vortex  cancels  the  dif¬ 
ferential  advection  by  the  vertical  shear  flow.  This  quasi¬ 
equilibrium  vortex  tilt  can  be  described  in  a  first  approxi¬ 
mation  by  the  balanced  dynamics  of  an  initially  barotropic, 
dry  TC-like  vortex  (Jones,  1995;  Reasor  et  al.,  2004).  In 
RMNV,x-  the  tilt  amplitude  decreases  while  maintaining  the 
quasi-equilibrium  tilt  direction.  The  TC  virtually  re-aligns 
within  approx.  24  h,  i.e.  after  this  time  the  tilt  amplitude  is 
comparable  to  the  tilt  amplitude  of  the  TC  in  the  no-shear 
benchmark  experiment  (see  Fig.  5  in  RMN10).  Once  the  tilt 
amplitude  has  decreased  considerably,  a  high  variability  in 
the  tilt  direction  ensues,  indicative  of  a  high-frequency  pre¬ 
cession  of  the  vortex  (Fig.  4a). 

Figures  3  and  4  demonstrate  that  the  general  characteris¬ 
tics  of  the  tilt  evolution  are  similar  in  all  of  our  experiments  in 
the  first  24-30  h  after  vertical  wind  shear  is  imposed.  All  TCs 
are  resilient  during  this  time  and  the  tilt  direction  is  predom¬ 
inantly  to  the  South  (180°).  In  CBLAST^,  the  maximum  tilt 
magnitude,  as  well  as  the  temporal  tilt  evolution,  is  very  sim¬ 
ilar  to  that  of  R.VIN(,8  •  The  only  exception  being  that  the  local 
maximum  at  23  h  is  more  pronounced  and,  thereafter,  the  tilt 
magnitude  remains  3-5  km  larger.  As  can  be  expected  from 
the  discussion  in  Sect.  3.1,  the  weaker  TCs  in  RMN54  and 
CBLAST54  clearly  exhibit  a  larger  tilt  magnitude  in  the  first 
24  h  than  the  more  intense  TCs  in  RMN68,  CBLAST(,8,  and 
ICE68-  Maximum  tilt  amplitudes  reach  40  km  in  RMN54  and 
30km  in  CBLAST54.  Both  TCs  eventually  re-align. 

The  most  striking  feature  in  the  tilt  evolution  is  the  dis¬ 
ruption  of  the  TC  vortex  after  30  h  in  ICEgs  (Fig.  3).  Until 
this  time  the  TC  tends  to  re-align  slowly:  the  local  tilt  maxi¬ 
mum  and  minimum  at  25  h  and  30  h,  respectively,  are  smaller 
than  their  respective  predecessors  at  9h  and  15  h.  At  30  h, 
intensity  has  decreased  by  over  30  m  s-1  as  compared  to  the 
time  when  vertical  wind  shear  is  imposed  and  has  reached  the 
smallest  absolute  value  in  all  experiments  (approx.  35  m  s_ , 
Fig.  2).  Based  on  the  discussion  in  Sect.  3.1,  it  is  plausible 
that  the  significantly  weakened  TC  can  no  longer  withstand 
the  differential  advection  associated  with  the  vertical  wind 
shear.  Note  that  the  distinct  weakening  in  the  first  20  h  pre¬ 
cedes  the  disruption  of  the  TC  by  more  than  10  h.  After  30  h, 
the  TC  remains  as  a  coherent  entity  with  smaller  vertical  ex¬ 
tent:  the  tilt  as  inferred  from  the  vorticity  centroid  at  9  km 
height  instead  of  10  km  height  does  not  exceed  35  km  (not 
shown).  Arguably,  the  survival  as  a  vertically  less  extensive 
system  allows  for  the  rapid  re-intensification  of  the  TC  after 
the  inflow  layer  0t  values  have  recovered  (Sect.  4.1). 

The  tilt  magnitude  in  ioICE<-,s  remains  below  10  km  and 
quickly  reduces  to  values  smaller  than  6  km  after  12  h.  De¬ 
spite  this  small  tilt  magnitude  a  southerly  tilt  direction  is  pre¬ 
ferred  until  late  in  the  experiment.  Although  the  evolution  of 
the  tilt  magnitude  is  similar  in  ioRMN^s,  the  tilt  direction 


in  10RMN68  exhibits  more  variability  after  30  h  when  the  tilt 
magnitude  is  small  (not  shown). 

3.4  SBC  formation  and  downward  flux  of  low-0e  air  into 
the  inflow  layer 

A  metric  for  the  downward  flux  of  low-0e  air  into  the  inflow 
layer  (DFX)  has  been  defined  in  RMN10  (their  Sect.  4.2.1) 
as14 

DF X  =  w-6'e,  (4) 

where  u>-  denotes  downward  vertical  motion  and  0'c  the  per¬ 
turbation  from  the  azimuthal  mean  of  0e.  The  quantity  DFX 
is  evaluated  at  1.5  km,  the  approximate  top  of  the  axisym- 
metric  inflow  layer. 

The  time-averaged  structure  of  the  SBC,  represented  by 
low-level  upward  motion,  and  the  distribution  of  DFX  shortly 
after  the  time  when  vertical  wind  shear  is  imposed  is  depicted 
in  Fig.  5.  In  this  early  stage  of  the  evolution,  the  general 
structure  of  the  SBC  and  the  distribution  of  DFX  is  consis¬ 
tent  in  all  experiments.  The  SBC  forms  in  the  downshear- 
right  quadrant15  extending  radially  outwards  to  150-200  km. 
Small  to  moderate  DFX  values  occur  over  a  broad,  banded 
region  downwind  of  the  SBC,  from  the  downshear  quadrant 
through  downshear-left,  and  into  the  upshear  quadrant.  In  all 
experiments,  azimuthal  vertical  cross-sections  (as  Fig.  18  in 
RMN10,  not  shown  here)  indicate  that  the  associated  down- 
drafts  form  underneath  the  helical  updrafts  of  the  SBC. 

It  is  of  interest  to  note  that  the  SBCs  in  RMN54,  ICE<-,8, 
and  10ICE68,  respectively,  extend  azimuthally  further  into 
the  upshear  quadrant  than  their  counterparts  in  RMN68, 
CBLAST68,  and  CBLAST54.  In  the  former  experiments,  the 
tilt  direction  at  early  times  exhibits  a  greater  downshear 
(westward)  component  than  in  the  latter  experiments  (Fig.  4). 
Thus,  the  low-level  vorticity  anomaly  associated  with  the 
tilted  vortex  can  be  expected  to  extend  into  the  upshear  quad¬ 
rant.  RMN10  have  shown  evidence  that  frictional  conver¬ 
gence  due  to  the  low-level  vorticity  anomaly  associated  with 
the  tilted  vortex  contributes  to  the  forcing  of  the  SBC.  The 
correlation  of  the  location  of  SBC  formation  with  the  tilt  di¬ 
rection  in  our  experiments  presented  here  provides  further 
support  for  RMNlO’s  interpretation. 

14Our  downward-flux  metric  DFX  is  virtually  identical  to  that 
part  of  the  flux  term  0'tw'  for  which  w'  is  negative  (see  RMN10, 
their  Sect.  4.2.1,  for  details). 

15  To  describe  the  azimuthal  location  of  relevant  features  in  this 
study  the  following  convention  is  used:  First,  the  centerline  of  the 
azimuthal  location  is  given  with  respect  to  the  shear  vector.  Then, 
the  azimuthal  extent  is  given,  e.g.  quadrant  or  semicircle.  For  ex¬ 
ample,  let  0°  denote  north  and  let  the  shear  vector  point  to  the  west 
(270°).  Then,  “downshear-right”  is  at  0°  and  the  “downshear-right 
quadrant”  denotes  the  region  between  315°  and  45°.  Note  that  this 
convention  differs  somewhat  from  the  one  that  is  frequently  used  in 
the  TC  community. 
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Fig.  5.  Downward  flux  of  8 'e,  DFX  (colour,  in  0.1  Kms  1 ,  see  text  for  definition),  at  the  top  of  the  axisymmetric  inflow  layer,  defined  as 
1 .5  km,  and  low-level  upward  motion  (thin  contour:  0.25  ms-1,  thick  contour:  1ms-1,  averaged  from  1 .25-2  km  height)  both  averaged  from 
l-6h  after  vertical  shear  is  imposed  for  (a)  RMNgg,  (b)  CBLASTgg,  (c)  RMN54,  (d)  CBLAST54,  (e)  ICEgg,  and  (f)  jglCEgg.  The  arrow 
indicates  the  direction  of  the  deep-layer  shear.  The  horizontal  scale  is  in  km. 


A  distinct  maximum  of  DFX  is  found  in  the  upshear  quad¬ 
rant  close  to  the  eyewall.  This  maximum  is  most  pronounced 
in  the  ICE  experiments  (Fig.  5e,  f).  In  RMN10,  this  DFX 
maximum  was  not  given  special  attention.  A  pronounced 
maximum  in  the  upshear  quadrant  occured  in  the  “20mps” 
case  of  RMN10  (their  Fig.  lOd).  RMN10  further  noted  a  lo¬ 
calized  region  of  downdrafts  in  the  upshear  region  close  to 
the  eyewall  after  the  TCs  in  the  “lOmps”  and  “15mps”  cases 
have  re-aligned  (their  Fig.  11a  and  associated  discussion). 

The  DFX  maximum  is  located  at  the  downwind  end  of  the 
banded  region  of  positive  DFX  values.  A  pronounced  6t  de¬ 
pression  is  found  at  the  top  of  the  inflow  layer  in  this  region 
(not  shown).  Therefore,  large  negative  values  of  0'  contribute 
to  high  values  of  DFX  (Eq.  4).  However,  downward  motion 
into  the  inflow  layer  is  particularly  strong  in  this  region  also 
(not  shown).  A  preliminary  examination  of  the  distributions 
of  vertical  motion,  hydrometeors,  and  0e  suggests  that  the 
downdrafts  close  to  the  eyewall  may  form  by  evaporation  of 
precipitation  falling  out  of  the  tilted  eyewall.  To  us,  however, 
it  is  not  clear  whether  the  more  pronounced  localized  max¬ 
imum  in  the  ICE  experiments  is  solely  due  to  the  inclusion 
of  ice  microphysics  or  to  what  extent  the  differences  in  the 
radial  structure  of  the  TC  vortex  contribute.  Clearly,  more  re¬ 


search  is  needed  to  clarify  the  formation  mechanism  of  these 
downdrafts,  as  well  as  the  relative  importance  of  the  confined 
DFX  maximum  and  the  broad  region  of  moderate  DFX  for 
TC  intensity  evolution. 

3.5  Vortex-tilt  induced  dynamic  instability? 

A  non-linear,  Kelvin-Helmholtz  type  instability  in  the  inner 
core  of  a  tilted,  dry  TC-like  vortex  has  been  found  in  a  nu¬ 
merical  experiment  by  Reasor  et  al.  (2004,  their  Sect.  5a).  In 
their  experiment,  this  instability  contributed  to  enhanced  lo¬ 
calized,  three-dimensional  mixing.  Such  mixing  is  one  path¬ 
way  to  dissipate  kinetic  energy  in  the  TC’s  inner  core.  The 
role  of  this  mixing  for  the  intensity  evolution  of  TCs  in  ver¬ 
tical  wind  shear  is  hitherto  unexplored. 

The  necessary  condition  for  such  an  instability  is  that  the 
gradient  Richardson  number,  7Z ,  drops  below  unity  (Abar- 
banel  et  al.,  1986).  Here,  TZ  is  calculated  as 


U  =  N2/«du/dz)2  +  (dv/dz)2),  (5) 
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Fig.  6.  Snapshots  of  the  gradient  Richardson  number  at  1625  m  height  and  3  h  after  vertical  shear  is  imposed  (shaded  for  values  below  1)  in 
(a)  RMNgg  and  (b)  ICEgg .  Upward  motion  is  depicted  in  black  contours  (thin  for  0.25  m  s-  1 ,  thick  for  1ms-1).  The  cyan  contour  denotes 
regions  of  strong  downward  flux  of  low-#e  air  into  the  inflow  layer  (DFX  at  1 .5  km  height,  the  contour  depicts  the  4.5  K  m  s- 1  isoline).  The 
horizontal  scale  is  in  km. 


where  N  is  the  (dry)  vertical  static  stability,  and  u  and  v  are 
the  radial  and  tangential  wind  components16,  respectively. 
The  tilt  of  the  vortex  may  lead  to  a  high  local  vertical  shear 
of  the  swirling  winds  in  the  inner  core.  Thus,  in  such  regions, 
the  gradient  Richardson  number  may  reach  small  values. 

In  all  of  our  (moist)  TC  experiments,  the  potential  for 
a  Kelvin-Helmholtz  type  instability  and  the  associated  three- 
dimensional,  turbulent  mixing  in  the  inner  core  is  an  ele¬ 
ment  of  vertical-shear  interaction.  Underneath  the  eyewall 
and  the  SBC,  TZ  <  1  is  an  ubiquitous  feature  within  the  layer 
of  strong,  frictionally-driven  inflow  (not  shown).  Above  this 
inflow  layer,  localized  regions  of  TZ  <  1  occur,  at  least  tem¬ 
porarily,  in  all  experiments  presented  in  this  study  within  the 
first  6h  after  shear  is  imposed.  These  regions  are  located 
within  the  eyewall  in  the  downshear-left  to  upshear  quad¬ 
rant  (exemplified  in  Fig.  6  for  RMNV,g  and  ICE(,g).  This  fea¬ 
ture  is  most  pronounced  in  RMNgg  and  least  pronounced 
in  RMN54,  as  might  be  expected  for  the  weaker  vortex  in 
RMN54.  The  regions  of  TZ  <  1  are  frequently  located  radi¬ 
ally  inwards  (within  20-30  km)  of  regions  of  high  DFX  (ex¬ 
emplified  in  Fig.  6).  Due  to  this  proximity,  it  seems  justi¬ 
fied  to  hypothesize  that  enhanced  mixing  associated  with  the 
tilt-induced  instability  locally  modifies  the  downward  flux  of 
lo\v-0e  air  into  the  inflow  layer.  A  quantitative  diagnostic  of 
this  mixing  process  is  deferred  to  future  work. 

3.6  Inflow  layer  0e  depression 

3.6.1  0e  depression  outside  of  the  eyewall 

A  notable  depression  of  0e  values  in  the  inflow  layer  outside 
of  the  eyewall  occurs  in  all  experiments  2-3  h  after  shear  is 
imposed.  A  snapshot  of  the  height-averaged  inflow  layer  9e 
distribution  representative  of  the  early  part  of  the  shear  ex- 

16Normally,  the  definition  of  TZ  involves  the  zonal  and  merid¬ 
ional  wind  components.  Our  calculation  of  TZ  is  equivalent  to  this 
definition. 


periments  (at  5  h)  is  presented  in  Fig.  7.  As  can  be  expected 
from  the  DFX  patterns  discussed  above,  the  inflow  layer  9e 
distribution  is  similar  in  all  experiments  at  this  time:  Low- 
9e  air  wraps  from  the  downshear  quadrant  into  the  upshear 
quadrant  cyclonically  inward17.  Also  evident  in  all  cases  is 
the  pronounced  wave-number  1  signature  of  the  SBC  in  the 
vertical  motion  field  at  low  levels  above  the  inflow  layer. 

To  facilitate  the  intercomparison  of  the  experiments  we 
consider  now  the  deviation  of  the  horizontal  distribution  of 
9e  relative  to  its  initial  azimuthal  average  at  the  time  when 
shear  is  imposed  in  the  respective  experiment  (Fig.  8).  This 
figure  depicts  time-averaged  values  from  4-9  h  after  shear  is 
imposed.  In  all  cases,  the  inflow  layer  9e  values  outside  of  the 
inner-core  updrafts  are  significantly  depressed.  A  depression 
of  up  to  15  K  can  be  found  in  RMN'(,g,  RMN54,  and  ICEgg. 
Minimum  9e  values  are  generally  found  to  the  southeast  of 
the  centre  (upshear  to  downshear-left).  This  location  is  con¬ 
sistent  with  the  DFX  maximum  in  this  region  and  the  contin¬ 
uous  depression  of  9e  values  along  the  spiralling  inflow  leg, 
as  hypothesised  by  RMN10  (their  Fig.  11a).  The  inflow  layer 
9e  depression  extends  into  the  downshear-right  quadrant  also, 
outwards  to  a  radius  of  100-150km.  In  the  absence  of  down¬ 
ward  flux  of  low  9e  in  this  region  (cf.  Fig.  5)  the  inflow  layer 
9e  values  gradually  recover  while  air  parcels  circle  from  the 
upshear  quadrant  to  the  upshear-right  quadrant. 

The  inflow  layer  9e  depression  outside  of  the  eyewall  in 
the  CBLAST  experiments  is  less  pronounced  than  in  the 
comparable  RMN  experiments.  While  the  minima  in  RMN 
exceed  —  13K,  the  minima  in  CBLAST  have  values  of  ap¬ 
prox.  — 10  K.  This  less  pronounced  9e  depression  is  consis¬ 
tent  with  smaller  average  DFX  values  in  CBLAS  l/,g  and 

^Approximately,  this  9e  depression  is  associated  with  a  depres¬ 
sion  in  the  water  vapor  mixing  ratio  by  4-5  g  kg- 1  and  a  depres¬ 
sion  in  the  potential  temperature  by  2-3  K.  The  spatial  distribution 
of  these  fields  for  ICEgg  and  RMNgg  at  the  time  shown  in  Fig.  7 
can  be  found  on  the  ACPD  webpage  in  the  authors’  response  to  the 
anonymous  reviewer  (item  3). 
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Fig.  7.  As  Fig.  5,  but  for  a  snapshot  of  inflow  layer  9e  (colour,  averaged  over  lowest  1  km)  and  low-level  updrafts  at  5  h  after  vertical  shear 
is  imposed  for  (a)  RMNgg,  (b)  CBLASTgg,  (c)  RMN54,  (d)  CBLAST54,  (e)  ICEgg,  and  (f)  iglCEgg.  The  depicted  time  is  representative  for 
the  early  part  of  the  experiments. 


CBLAST54  as  compared  to  RMN^g  and  RMN54,  respec¬ 
tively  (Fig.  5).  The  smaller  DFX  values  in  CBLAST  are  asso¬ 
ciated  with  a  weaker  general  downdraft  activity  (see  Sect.  3.7 
and  Fig.  10  for  supporting  evidence). 

In  the  ICE  experiments,  the  spatial  pattern  of  the  9e  de¬ 
pression  is  radially  more  confined  than  in  the  warm-rain  ex¬ 
periments.  This  is  consistent  with  the  DFX  pattern.  As  dis¬ 
cussed  at  the  end  of  Sect.  3.4,  it  is  not  evident  to  us  to  what 
extent  this  greater  radial  confinement  can  be  attributed  to  dif¬ 
ferences  in  the  cloud  microphysics  or  to  differences  in  the 
radial  structure  of  the  TC  vortices. 

3.6.2  0e  depression  underneath  the  inner-core  updrafts 

An  axisymmetric  view  of  the  inflow  layer  and  eyewall  9e  dis¬ 
tribution  is  presented  in  RMN10.  In  this  axisymmetric  view 
a  tongue  of  low-0e  air  effectively  intrudes  into  the  inflow 
layer  from  above  and  radially  inwards  towards  the  eyewall. 
The  9e  values  of  air  rising  in  the  eyewall  of  the  sheared  TCs 
are  several  degrees  smaller  than  in  the  experiment  without 
vertical  shear  (Fig.  9  in  RMN10).  The  same  general  charac¬ 
teristics  of  the  axisymmetric  inflow  layer  and  adjacent  eye- 
wall  structure  are  found  in  the  additional  experiments  consid¬ 
ered  in  this  study  (not  shown).  When  viewed  from  the  sim¬ 


plified  perspective  of  an  idealised  Carnot-cycle,  a  weakening 
of  the  sheared  TCs  can  be  expected.  As  noted  in  RMN10, 
however,  it  is  an  open  question  whether  an  axisymmetric 
viewpoint  is  sufficient  in  the  presence  of  pronounced  wave- 
number  1  asymmetries  in  the  inflow  layer  9e  distribution  and 
the  eyewall  updrafts  of  a  TC  in  vertical  wind  shear. 

In  the  present  study  we  complement  the  axisymmetric 
view  by  the  examination  of  the  horizontal  distribution  of  in¬ 
flow  layer  9t  and  low-level  eyewall  updrafts.  For  the  follow¬ 
ing  discussion  the  eyewall  updrafts  are  defined  as  the  con¬ 
tiguous  region  within  the  1ms-1  vertical  velocity  isotach 
near  the  centre  in  Fig.  8.  We  consider  the  9t  values  within 
this  region  to  be  a  good  approximation  for  the  9e  values  of  air 
rising  out  of  the  boundary/  sub-cloud  layer  into  the  eyewall 
updrafts.  Upon  examining  the  results,  it  is  obvious  that  in 
all  cases,  the  inner-core  updrafts  exhibit  a  pronounced  wave- 
number  1  asymmetry  and  that  the  eyewall  is  not  closed. 

By  visual  inspection  we  estimate  the  average  9e  depres¬ 
sion  underneath  the  eyewall  (averaged  from  4—9  h)  to  be  0  K 
for  CBLAST54,  2-3  K  for  RMN68  and  CBLAST68,  3-5  K  for 
RMN54  and  10ICE68,  and  6-7  K  for  ICE(,g.  These  values  cor¬ 
relate  well  with  the  intensity  evolution  of  the  respective  TCs. 
The  TC  in  CBLAST54  does  not  weaken  until  7  h  after  shear 
is  imposed  (Fig.  2c).  The  TCs  in  RMNe8  and  CBLAST68 
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Fig.  8.  Depression  of  inflow  layer  9e  (averaged  over  the  lowest  1  km  and  from  4-9  h)  relative  to  the  axisymmetric  mean  at  the  time  when 
vertical  wind  shear  is  imposed  (colours).  Thin  contours  depict  the  radial  flow  averaged  over  the  lowest  1  km  (every  4ms-1,  starting  at 
2ms-1,  dashed  for  negative  values  (outflow)).  Thick  contour  denotes  low-level  updrafts  (lms-  ,  averaged  from  1.25-2 km  height),  (a) 
RMN68,  (b)  CBLAST68,  (c)  RMN54,  (d)  CBLAST54,  (e)  ICE68,  and  (f)  10ICE68. 


show  a  very  similar  weakening  in  the  first  12  h  while  RMN54 
and  i0ICE68  exhibit  slightly  faster  weakening  rates.  The  pro¬ 
nounced  and  rapid  weakening  of  the  TC  in  ICE()8  is  associ¬ 
ated  with  the  most  pronounced  depression  of  9e  underneath 
the  eyewall.  This  correlation  supports  RMNlO’s  result  that 
the  intensity  evolution  of  resilient  TCs  in  vertical  shear  is 
governed  to  first  order  by  the  systematic  depression  of  9t  in 
the  inflow  layer. 

The  inflow  layer  9e  depression  outside  of  the  eyewall  in 
CBLAST()8  is  less  pronounced  than  in  RMN^x-  However,  the 
9e  depression  underneath  the  eyewall  is  very  similar  in  both 
experiments.  Consistently,  the  intensity  decrease  is  very  sim¬ 
ilar  in  CBLASTfix  and  RMN  6  8  in  the  first  12  h.  Therefore, 
it  can  be  inferred  that  the  replenishment  of  inflow  air  spi¬ 
ralling  towards  the  eyewall  is  more  complete  in  RMN  than  in 
CBLAST.  As  discussed  in  the  Introduction,  the  smaller  value 
of  Cd  in  CBLAST  may  lead  to  weaker  frictional  inflow  as 
compared  to  RMN  and  thus  to  a  longer  replenishment  time. 
Here,  this  effect  of  a  smaller  Co  value  is  apparently  over¬ 
ruled  by  the  reduced  extraction  of  heat  from  the  ocean  sur¬ 
face  due  to  a  smaller  value  of  Ck ■  It  is  difficult  to  examine 
this  aspect  in  RMN54  and  CBLAST54  because  the  9e  depres¬ 


sion  underneath  the  eyewall  is  different  at  this  early  time  in 
the  experiment. 

The  replenishment  time  of  air  with  depressed  9e  val¬ 
ues  may  depend  considerably  on  the  azimuthal  location. 
The  radial  velocity  averaged  over  the  lowest  1  km  exhibits 
a  pronounced  wave-number  1  asymmetry  in  all  experiments 
(Fig.  8).  In  the  downshear  semicircle,  the  maximum  height- 
averaged  radial  inflow  is  approx.  20  ms-1  while  in  the  up- 
shear  semicircle  the  height-averaged  radial  velocity  is  very 
small  or  is  even  directed  outwards  (e.g.  in  RMN54  (Fig.  8c)). 
Figure  8  indicates  that  the  bulk  of  the  depressed  inflow  layer 
9e  spirals  around  the  centre  before  entering  the  inner  core 
updrafts  in  the  downshear-right  quadrant  (e.g.  Fig.  8a,  f).  In 
contrast,  in  RMN54  and  ICE^s,  an  additional,  much  more  di¬ 
rect  pathway  for  low-()c  air  from  the  downshear  quadrant  to 
enter  the  eyewall  updrafts  in  the  downshear  to  downshear-left 
quadrant  is  indicated  (Fig.  8c,  e).  The  replenishment  of  in¬ 
flow  layer  9e  along  this  more  direct  pathway  can  be  expected 
to  be  less  complete  than  along  the  longer  inflow  spiral.  In 
the  presence  of  a  pronounced  wave-number  1  asymmetry  of 
the  inflow,  it  stands  to  reason  that  the  (azimuthal)  location 
of  the  flushing  of  the  inflow  layer  with  low-/9e  air  may  be  of 
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Fig.  9.  Contoured  frequency  by  altitude  diagrams  of  vertical  motion  within  160  km  radius  of  the  centre.  Vertical  motion  is  binned  into 
intervals  of  0.5  ms~  ,  centred  on  the  values  given  on  the  r-axis.  The  colour-coded  value  denotes  the  logarithm  of  the  frequency.  With  each 
contour  the  frequency  decreases  by  a  factor  of  1/e.  Before  taking  the  logarithm  of  the  frequency  a  small  positive  value  (10-5  ss  e~ 1 15)  has 
been  added  to  ensure  that  the  logarithm  is  well  defined  everywhere.  For  the  no-shear  cases  of  RMN  (a)  and  ICE  (b)  the  6  h  time  period  after 
vertical  wind  shear  is  imposed  in  the  respective  shear  runs  is  considered.  For  RMNgg  (c)  and  ICEgg  (d)  the  times  4-9  h  are  considered. 


importance  for  the  intensity  evolution  also,  and  not  only  the 
magnitude  of  the  0e  depression  in  the  azimuthal  mean. 

3.7  Vertical  motion  characteristics 

The  formation  of  the  SBC  and  the  associated  downdrafts 
that  flush  low-tie  air  into  the  inflow  layer  are  two  central  as¬ 
pects  of  RMNlO’s  framework.  Here,  we  examine  the  gen¬ 
eral  characteristics  of  vertical  motion  in  our  suite  of  exper¬ 
iments  using  so-called  contoured  frequency  by  altitude  di¬ 
agrams  (CFADs,  Yuter  and  Houze,  1995).  These  diagrams 
depict  the  frequency  of  the  occurrence  of  vertical  motion  of 
a  given  magnitude,  here  binned  in  0.5  m  s~ 1  intervals,  within 
a  given  region  and  at  a  given  altitude.  It  is  very  difficult  to 
focus  specifically  on  the  SBC  and  on  regions  where  shear- 
induced  downdrafts  form.  Therefore,  we  consider  all  vertical 
motion  within  a  radius  of  160  km. 

A  comparison  of  the  RMN  and  ICE  experiments  is  pre¬ 
sented  in  Fig.  9.  Without  vertical  wind  shear  (Fig.  9a,  b), 
strong  up-  and  downdrafts  occur  more  frequently  in  the  ICE 
experiment.  Under  the  influence  of  vertical  shear  (Fig.  9c, 
d),  the  frequency  distribution  broadens  considerably  below 
5-6  km  as  compared  to  the  respective  no-shear  experiment. 
For  both  the  warm  rain  and  ice  microphysics  configurations, 
vertical  wind  shear  increases  the  frequency  of  the  strongest 
up-  and  downdrafts  at  midlevels  and  below.  At  upper  lev¬ 


els  the  increase  of  downdraft  frequency  is  less  pronounced 
and  the  frequency  of  strong  updrafts  tends  to  decrease,  ex¬ 
cept  for  the  strongest  updrafts  in  ICEgs-  It  is  evident  also 
that  strong  downdrafts  below  5  km  occur  considerably  more 
frequently  in  ICEgg  than  in  RMNgs  (cf.  Fig.  9c,  d).  The 
more  pronounced  downdraft  activity  in  the  ICE  experiments 
is  consistent  with  the  tenet  that  the  cooling  associated  with 
the  melting  of  ice  particles  may  initiate  or  strengthen  exist¬ 
ing  downdrafts  (Srivastava,  1987). 

Whithout  imposed  vertical  shear,  the  results  in  CBLAST 
are  very  similar  to  those  in  RMN  (not  shown).  The  im¬ 
posed  vertical  wind  shear  leads  to  a  broadening  of  the  fre¬ 
quency  distribution  in  CBLAST  also  (not  shown).  The  dif¬ 
ferences  between  the  CBLAST  and  RMN  shear  experiments 
are  less  pronounced  than  between  ICE  and  RMN.  The  dif¬ 
ferences  are  depicted  in  Fig.  10  as  the  ratio  of  the  fre¬ 
quencies  in  RMN68  and  CBLAST68  (Fig.  10a)  and  RMN54 
and  CBLAST54  (Fig.  10b).  Apparently,  the  frequency  of  the 
strongest  downdrafts  is  more  than  275  %  greater  in  RMN 
than  in  CBLAST.  Strong  updrafts  tend  to  be  more  frequent 
in  RMN  also,  in  particular  so  in  RMN54  as  compared  to 
CBLAST54. 

Differences  in  the  thermodynamic  properties  of  the  am¬ 
bient  air  in  which  downdrafts  occur  may  contribute  to  the 
observed  differences  in  downdraft  activity.  The  strength  of 
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Fig.  10.  As  Fig.  9,  but  for  the  ratio  of  vertical  motion  frequency  of  RMNgg/CBLASTgg  (a)  and  RMN54/CBLAST54  (b). 


a  downdraft  depends,  inter  alia,  on  the  buoyancy  difference 
between  the  descending  air  parcel  and  the  ambient  air  along 
the  parcel’s  downward  path.  A  vertically  integrated  mea¬ 
sure  of  this  buoyancy  difference  is  the  downdraft  convec¬ 
tive  available  potential  energy  (DCAPE,  Emanuel,  1994,  his 
Sect.  6.3).  The  derivation  of  DCAPE  assumes  the  following 
two-stage  process.  First,  evaporation  of  precipitation  cools  an 
air  parcel  isobarically  until  saturation  is  reached.  Then,  the 
parcel  descends  pseudo-adiabatically  with  just  enough  evap¬ 
oration  to  keep  the  parcel  saturated. 

We  have  calculated  DCAPE  using  the  azimuthal  mean 
soundings  at  the  time  when  shear  is  imposed.  Because  it  is 
not  entirely  clear  which  radius  is  most  appropriate  to  ex¬ 
amine  the  downdrafts  underneath  the  SBC  we  considered 
a  number  of  different  radii  and  radial  averages  between 
60km  and  150km  (not  shown).  For  none  of  these  regions, 
could  we  find  a  clear  relationship  between  the  DCAPE  val¬ 
ues  and  the  downdraft  activity,  or  between  the  DCAPE  values 
and  the  downward  flux  of  low-0e  air  into  the  inflow  layer.  In 
particular,  DCAPE  values  below  500  hPa  in  CBLAST54  are 
higher  than  in  RMN54,  in  contrast  to  the  more  pronounced 
downdraft  activity  in  RMN54. 

Apparently,  consideration  of  solely  thermodynamic  prop¬ 
erties  of  the  air  in  which  the  downdrafts  occur  is  insufficient 
to  explain  the  observed  differences  in  downdraft  activity.  We 
propose  the  following  alternative  explanation  that  focuses  on 
the  dynamic  coupling  of  the  convection  with  the  frictional 
boundary  layer.  First,  we  note  that  the  weaker  downdrafts 
in  CBLAST68  and  CBLAST54  are  associated  with  weaker 
updraft  activity.  This  may  be  expected  because  a  prominent 
downdraft  formation  mechanism  is  the  evaporation  of  pre¬ 
cipitation  that  is  formed  in  the  updrafts  and  falls  into  unsatu¬ 
rated  air.  Both  the  smaller  value  of  Cd  and  C k  in  CBLAST 
may  contribute  to  weaker  updrafts.  Ascent  out  of  the  bound¬ 
ary  layer  forced  by  frictional  convergence,  either  underneath 
the  eyewall  or  underneath  the  vorticity  anomaly  associated 
with  the  tilt  of  the  outer  vortex  (RMN10,  their  Fig.  17),  can 
be  expected  to  decrease  for  decreasing  Cd-  For  example,  in 
simple  Ekman  pumping  theory  the  vertical  motion  at  the  top 
of  the  boundary  layer  is  proportional  to  ~JCd  (e.g.,  Holton, 


2004,  his  Eq.  5.38).  The  smaller  value  of  Ck  in  CBLAST 
decreases  the  rate  by  which  heat  is  extracted  from  the  ocean 
surface.  These  reduced  surface  fluxes  destabilise  the  bound¬ 
ary  layer  in  CBLAST  to  a  lesser  extent  than  in  RMN.  Less 
vigorous  updrafts  can  be  expected  to  occur  in  such  a  more 
stable  stratification.  Our  foregoing  argument  is  supported  by 
the  results  of  the  CBLAST  and  RMN  no-shear  experiments 
(not  shown).  Without  vertical  shear,  the  activity  of  transient 
rain  bands  and  the  associated  localized  downward  flux  of  low 
9e  into  the  inflow  layer  in  CBLAST  are  less  pronounced  than 
in  the  no-shear  RMN  experiment  also. 

In  all  experiments,  vertical  wind  shear  increases  the  fre¬ 
quency  of  moderate  to  strong  downdrafts.  It  is  clear  that 
CBLAST  has  the  lowest  frequency  of  strong  downdrafts  and 
that  ICE  has  the  highest  frequency.  The  downdraft  activity  as 
seen  in  the  CFADs  is  consistent  with  the  magnitude  of  the 
inflow  layer  9e  depression  outside  of  the  eyewall  (cf.  Fig.  8). 
Therefore,  the  examination  of  the  CFADs  provides  additional 
support  for  the  hypothesis  that  the  distinct  depression  of  the 
inflow  layer  9e  values  are  caused  by  downdrafts  induced  by 
the  imposed  vertical  wind  shear. 

4  A  closer  look  at  the  ICKciS  experiment 

The  weakening  of  the  TC  in  ICE68  is  much  more  pronounced 
than  in  the  comparable  warm-rain  experiments  RMN68  and 
CBLAST68-  Nevertheless,  in  the  previous  sections  it  has 
been  shown  that  the  general  structural  changes  in  the  first 
20  h  are  very  similar  in  all  experiments.  We  now  take  a  closer 
look  at  ICE(38  to  further  confirm  that  the  intensity  evolution 
in  this  case  is  consistent  with  the  shear-induced  thermody¬ 
namic  modification  of  the  inflow  layer.  First,  the  inflow  layer 
9e  distribution  is  examined  during  periods  of  distinct  inten¬ 
sity  changes.  Then,  we  will  provide  evidence  based  on  ax- 
isymmetric  spindown  theory  that  the  rapid  weakening  in  the 
first  20  h  is  consistent  with  the  spindown  of  the  TC  vortex  due 
to  divergence  above  the  inflow  layer  associated  with  a  signif¬ 
icant  decrease  of  the  inner-core  convective  mass  flux. 


Atmos.  Chem.  Phys.,  13,  327-346,  2013 


www.atmos-chem-phys.net/13/327/2013/ 


M.  Riemer  et  al.:  Thermodynamic  impact  of  vertical  wind  shear  on  the  inflow  layer 


343 


339  343  347  351  355  359  363  367  371  375  [K] 


339  343  347  351  355  359  363  367  371  375  [K] 


339  343  347  351  355  359  363  367  371  375  [K] 


Fig.  11.  As  Fig.  7,  but  for  ICE68  at  (a)  7  h,  (b)  20  h,  (c)  23  h,  (d)  40  h,  (e)  45  h,  and  (f)  64  h. 


4.1  Intensity  evolution  and  inflow  layer  6e 

It  has  been  shown  in  Sect.  3.6  (Fig.  7e)  that  the  rapid  weaken¬ 
ing  is  associated  with  a  significant  depression  of  inflow  layer 
0e  soon  after  shear  is  imposed.  Here  we  emphasize  that  a  sig¬ 
nificant  part  of  the  depressed-(9c  air  apparently  intrudes  from 
the  downshear  quadrant  into  the  inner-core  updrafts  on  a  di¬ 
rect  pathway  (Fig.  11a,  at  7  h).  Over  half  of  the  area  under¬ 
neath  the  inner-core  updrafts  at  this  time  has  0e  values  lower 
than  355  K,  as  compared  to  approx.  363  K  before  shear  is  im¬ 
posed  (not  shown).  Along  with  this  significant  frustration  of 
the  0e  values  in  the  inner-core  region  the  vertical  mass  flux 
at  midlevels  (5  km)  decreases  in  the  first  10  h  by  60%  from 
50  x  106  kg  s-1  to  approx.  20  x  106  kg  s-1  (not  shown). 

We  consider  now  the  two  distinct,  short-lived,  local  inten¬ 
sity  maxima  at  22  h  and  44  h,  prominent  features  of  the  in¬ 
tensity  evolution  in  ICE68  (Fig.  2b,  the  maxima  are  found 
at  64  h  and  86  h  because  this  figure  includes  the  42  h  spinup 
period).  Both  periods  of  temporary  re-intensification  are  pre¬ 
ceded  by  relatively  high  inflow  layer  6e  values  in  the  downs¬ 
hear  to  downshear-right  quadrants  while  the  very  pronounced 
0e  depression  at  these  times  is  mostly  confined  to  the  up- 
shear  semicircle  (at  20 h  (Fig.  lib)  and  at  40 h  (Fig.  lid)). 
Due  to  the  distinct  wave -number  1  inflow  asymmetry  dis¬ 
cussed  previously  (Fig.  8),  it  can  be  expected  that,  at  these 
times,  the  eyewall  updrafts  are  fed  to  a  large  extent  by  the 


relatively  high-0e  air  in  the  downshear  semicircle.  Eyewall 
0e  values  may  then  temporarily  increase  and  the  TC  may  in¬ 
tensify.  The  temporary  intensification  ensues  until  the  very 
low  0e  air  from  the  upshear  semicircle  wraps  around  the  cen¬ 
tre  and  limits  considerably  the  supply  of  relatively  high  0e 
(exemplified  at  23  h  (Fig.  11c)  and  at  45  h  (Fig.  1  le)). 

At  the  end  of  our  experiment  a  robust  re-intensification  is 
indicated  in  the  intensity  time  series  (Fig.  2b).  As  in  the  ex¬ 
periments  examined  in  RMN10,  this  final  re-intensification 
is  associated  with  an  almost  complete  recovery  of  the  inflow 
layer 0e  (Fig.  Ilf). 

As  indicated  in  Fig.  llb-e,  low-0e  air  in  the  inflow  layer 
may  wrap  around  the  centre  into  the  downshear-right  quad¬ 
rant.  This  low- fic  air  tends  to  suppress  the  convection  of  the 
SBC.  Furthermore,  after  20  h  in  ICE()8,  the  vortex  tilt  is  no 
longer  in  the  approximate  left-of-shear  equilibrium  direction, 
but  starts  to  precess.  In  accordance  with  this  precession,  the 
vorticity  anomaly  associated  with  the  tilt  of  the  outer  vortex 
then  starts  to  move  also  and  thus  no  longer  provides  persis¬ 
tent  forcing  for  the  SBC  in  the  downshear-right  quadrant.  It 
is  interesting  to  note  that  a  narrow  band  of  convection  prefer¬ 
entially  occurs  in  the  downshear  quadrant  at  these  later  times 
(Fig.  1  lb-e).  Frequently,  one  observes  significant  downward 
flux  of  low  0e  associated  with  this  remnant  SBC  (not  shown). 
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4.2  A  timescale  for  vortex  spindown 

Within  the  simplified  framework  of  an  idealised  Carnot  cy¬ 
cle,  and  based  on  the  idea  that  environmental  vertical  wind 
shear  leads  to  the  frustration  of  such  a  thermodynamic  cy¬ 
cle,  Tang  and  Emanuel  (2010)  have  derived  a  theory  for  the 
storm’s  reduced  steady-state  intensity.  Such  a  theory,  how¬ 
ever,  does  not  provide  information  about  the  temporal  evo¬ 
lution  of  storm  intensity,  i.e.  how  rapidly  a  TC  weakens. 
A  well-known  concept  for  the  weakening  of  a  TC  is  the  ax- 
isymmetric,  frictionally-induced  spindown  of  the  TC  vortex. 
Frictional  convergence  forces  ascent  out  of  the  inflow  layer. 
When  the  vertical  convective  mass  flux  is  smaller  than  the 
frictionally  forced  vertical  mass  flux,  some  of  the  air  has 
to  diverge  above  the  inflow  layer  due  to  mass  continuity 
(Willoughby,  1979).  Under  the  assumption  that  air  parcels 
approximately  conserve  their  angular  momentum  above  the 
frictional  inflow  layer,  this  divergence  implies  a  decay  of  the 
tangential  winds,  i.e.  the  TC  weakens. 

A  timescale  for  the  frictional  spindown  of  an  axisymmet- 
ric  vortex  without  convective  mass  flux  has  been  derived  by 
Elliassen  and  Lystad  (1977).  Their  model  assumes  a  neutral 
stratification  above  the  layer  of  strong  frictional  inflow  which 
is  a  reasonable  assumption  in  the  inner  core  of  a  TC  because 
isosurfaces  of  9e  and  angular  momentum  are  approximately 
congruent  in  this  region.  Montgomery  et  al.  (2001)  have  ver¬ 
ified  by  numerical  experiments  that  the  theory  gives  reason¬ 
able  results  for  the  spindown  of  TC-like  vortices,  even  be¬ 
yond  the  formal  validity  of  the  theory.  The  half-time  /half, 
i.e.  the  time  that  is  required  to  reduce  the  initial  tangential 
velocity  v  by  a  factor  of  2,  is  given  by  Elliassen  and  Lystad’ s 
theory  as 

H  —  h 

fiialf =  yy  ■  (6) 

X  CDv 

Here,  H  is  the  fluid  depth,  h  the  height  of  the  boundary  layer, 
and  x  a  reduction  factor  defined  by  Elliassen  and  Lystad  as 
the  ratio  of  the  tangential  winds  at  the  top  of  the  surface  layer 
to  the  tangential  winds  at  the  top  of  the  Ekman  layer.  Esti¬ 
mating  H  —  h—  104  m,  x  —  0-8,  and  using  Cd  —  3.8  x  10-3 
from  Deacon’s  formula  with  r  =  68ms_1  yields  a  half-time 
/fiai  r  —  16.8  h. 

The  weakening  of  the  TC  in  ICEgg  is  remarkable:  intensity 
decreases  from  68  m  s_1  at  42  h  to  33  m  s-1  at  68  h  (Fig.  2b), 
corresponding  to  a  weakening  rate  of  35 ms-1  in  26 h;  or 
a  half-time  of  approx.  25  h.  This  value  is  consistent  with 
a  partial  shutdown  of  the  convective  mass  flux.  As  noted 
above,  concurrently  with  the  frustration  of  the  inflow  layer 
9e  and  the  rapid  weakening,  the  vertical  mass  flux  at  5  km 
height  decreases  by  approx.  60%.  Consequently,  it  can  be 
argued  that  the  significant  reduction  of  the  convective  verti¬ 
cal  mass  flux  by  the  severe  depression  of  the  inflow  layer  0e 
values  broadly  governs  the  rapid  weakening  in  ICE^s- 


5  Conclusions 

In  the  light  of  a  new  framework  developed  in  a  recent  study 
(RMN10),  the  current  study  examines  further  the  evolution 
of  TCs  exposed  to  environmental  vertical  wind  shear.  The 
new  framework  has  been  developed  on  the  basis  of  idealized 
numerical  experiments  for  which  we  have  deliberately  cho¬ 
sen  a  simple  set  of  physical  parameterisations  to  minimize 
the  complexity  of  the  problem.  In  the  current  study,  we  as¬ 
sess  the  robustness  of  our  previous  results  in  a  more  realistic 
and  representative  experimental  setup  by  surveying  and  di¬ 
agnosing  five  additional  numerical  experiments. 

The  modifications  of  the  experimental  setup  focus  on  pa¬ 
rameters  that  are  important  for  structure  change  and  inten¬ 
sity  evolution  of  numerically  simulated  TCs  and  the  storms’ 
evolution  in  vertical  wind  shear,  discussed  in  some  detail  in 
the  Introduction  in  Sect.  1.2.  Specifically,  in  comparison  to 
RMN10,  we  here  consider  weaker,  but  still  mature  TCs;  em¬ 
ploy  values  for  the  exchange  coefficients  of  surface  heat  and 
momentum  fluxes  that  are  in  better  agreement  with  recent 
observational  results;  and  complement  our  previous  warm- 
rain  experiments  with  experiments  that  employ  an  ice  micro¬ 
physics  scheme. 

Several  environmental  factors,  such  as  the  vertical  profiles 
of  moist  entropy  and  the  environmental  winds,  likely  play  a 
role  for  the  evolution  of  TCs  in  vertical  wind  shear  also.  Fur¬ 
thermore,  the  motion  vector  of  the  TCs  in  all  of  our  experi¬ 
ments  is  approx.  30°-45°  left  of  the  shear  vector  (not  shown) 
while  TCs  in  the  real  atmosphere  generally  feature  a  more 
pronounced  and  more  variable  angle.  The  role  of  the  envi¬ 
ronmental  factors  and  of  storm  motion  constitutes  an  impor¬ 
tant  topic  for  future  research.  It  is  hoped  that  the  framework 
presented  in  RMN10,  Riemer  and  Montgomery  (2011)  and 
in  the  current  study  provides  a  helpful  conceptual  basis  for 
such  research. 

The  ability  of  TC-like  vortices  to  withstand  the  differential 
advection  by  vertical  wind  shear  and  to  eventually  re-align 
depends  on  the  radial  profile  of  the  swirling  winds  (discussed 
in  detail  in  Sect.  3.1).  Some  differences  in  this  radial  struc¬ 
ture  arise  in  the  individual  experiments  before  vertical  wind 
shear  is  imposed.  The  most  notable  differences  are  found  in 
the  experiment  with  ice  microphysics  in  which  the  TC,  inter 
alia,  features  a  radius  of  maximum  winds  that  is  approx.  25  % 
smaller  than  in  the  warm-rain  experiments.  Differences  in  the 
radial  vortex  structure  complicate  the  attribution  of  differ¬ 
ences  in  the  TC  evolution  to  the  inclusion  of  ice  microphysics 
alone. 

The  differences  in  the  vortex  structure  notwithstanding,  in 
all  experiments  the  TCs  exhibit  the  same  general  structural 
changes  when  exposed  to  vertical  wind  shear.  The  structural 
changes  are  consistent  with  RMNlO’s  new  framework  and 
comprise  a  vortex  tilt  predominantly  to  the  downshear-left 
for  20-24  h,  the  formation  of  a  pronounced  wave-number  1 
convective  asymmetry  outside  of  the  eyewall  (reminiscent  of 
a  stationary  band  complex),  and  downward  flux  of  low-0e  air 
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into  the  inflow  layer  in  the  downshear-left  semicircle.  A  sig¬ 
nificant  depression  of  the  inflow  layer  9e  values  (O  (10  K)) 
occurs  in  all  experiments. 

Furthermore,  for  all  experiments  the  potential  for  a  non¬ 
linear,  Kelvin-Helmholtz  type  instability  in  the  inner  core  of 
the  tilted  TC-vortex  is  diagnosed.  Such  an  instability  would 
contribute  to  enhanced  dissipation  of  kinetic  energy  and  may 
modify  the  flux  of  low-tie  air  into  the  inflow  layer.  However, 
a  quantitative  analysis  of  the  mixing  process  associated  with 
the  potential  instability  is  deferred  to  future  work. 

Two  additional  diagnostics  are  considered  in  the  current 
study  that  have  not  been  considered  in  RMN10,  namely  the 
9e  values  underneath  the  asymmetric  inner-core  updrafts,  and 
so-called  contoured  frequency  by  altitude  diagrams  to  exam¬ 
ine  the  general  characteristics  of  vertical  motion.  The  mag¬ 
nitude  of  the  9e  depression  underneath  the  asymmetric  inner- 
core  updrafts  early  in  the  respective  shear  experiment  corre¬ 
lates  well  with  the  ensuing  weakening  of  the  TC,  supporting 
the  idea  that  it  is  the  depression  of  the  inflow  layer  9e  that 
governs  the  intensity  evolution  to  first  order.  For  all  experi¬ 
ments,  vertical  wind  shear  increases  the  frequency  of  strong 
downdrafts  as  compared  to  the  no-shear  experiments.  The 
magnitude  of  the  9e  depression  outside  of  the  eyewall  is  con¬ 
sistent  with  the  frequency  of  strong  downdrafts,  supporting 
the  idea  that  it  is  the  shear-induced  downdrafts  that  flush  low- 
9e  air  into  the  inflow  layer.  Consequently,  both  additional  di¬ 
agnostics  corroborate  RMNlO’s  framework. 

In  one  of  our  experiments  employing  ice  microphysics, 
a  very  rapid  weakening  of  the  TC  occurs.  This  intensity 
evolution  is  distinct  from  the  intensity  evolution  in  our 
warm-rain  experiments.  Therefore,  a  question  may  natu¬ 
rally  arise  whether  ice  microphysics,  as  represented  in  the 
RAMS  model,  introduce  an  important  mechanism  for  inten¬ 
sity  change  that  is  not  considered  in  RMNlO’s  framework. 
As  noted  above,  the  structural  changes  caused  by  the  vertical 
wind  shear  in  the  experiment  with  ice  microphysics  are  very 
similar  to  those  in  the  warm-rain  experiments.  In  our  suite  of 
experiments,  the  general  downdraft  activity,  the  flux  of  low- 
9e  air  into  the  inflow  layer,  and  the  9e  depression  underneath 
the  inner-core  updrafts  are  all  the  most  pronounced  in  the 
ice  experiment.  The  inner-core  9e  depression  that  accompa¬ 
nies  the  initial  rapid  weakening  is  more  than  twice  as  high 
as  that  in  the  comparable  experiments  employing  warm-rain 
microphysics  (6-7  K  as  compared  to  2-3  K).  Axisymmetric 
spindown  theory  (Elliassen  and  Lystad,  1977;  Montgomery 
et  al.,  2001)  further  supports  the  hypothesis  that  the  rapid 
initial  weakening  is  due  to  a  substantial  decrease  of  the  con¬ 
vective  mass  flux  associated  with  the  significant  depression 
of  the  inflow  layer  9t  by  shear-induced  downdrafts.  In  addi¬ 
tion,  we  have  documented  a  close  relationship  between  the 
evolution  of  the  inflow  layer  9e  values  and  two  transient  in¬ 
tensity  maxima  that  occur  at  22  h  and  44  h  after  vertical  wind 
shear  is  imposed.  The  final  re-intensification  of  the  TC  in  this 
experiment,  as  in  RMN10,  is  associated  with  a  near-complete 
replenishment  of  the  inflow  layer  9e  values  also. 


Based  on  the  evidence  presented  here,  we  conclude  that 
RMNlO’s  framework  for  the  intensity  modification  of  TCs 
in  vertical  wind  shear  is  robust  to  the  presented  changes  in 
the  setup  of  the  numerical  experiment.  We  advocate  further 
study  and  tests  of  this  framework  with  case  studies  and  real 
atmospheric  data. 
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